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Research initiated under contract N00014-78-C-0889 is most conveniently 
summarized in the form of four appendices. First, in Appendix I, the five 
different areas of research on water waves are briefly summarized. In two of 
these areas (evolution of long internal waves and evolution of envelope solitons) 
work has not been completed; some results to date are described. In the remaining 
three areas, final manuscripts have been completed and are presented in Appendices 
II, III, and IV. In particular, Appendix II is a copy of a paper entitled 
'Baroclinic Tsunami Generation" to appear in the Journal of Physical Oceanography, 
September 1980; Appendix III is a copy of a paper entitled "Long Waves Generated 
by Complex Bottom Motions" to appear in Proceedings, 17th Conference on Coastal 
Engineering; and Appendix IV is a copy of lecture notes on small-scale ocean 
waves presented at the summer (1980) course "Topics In Ocean Physics" held at 
the International School of Physics "Enrico Fermi", Varenna, Italy. 
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Research Summary 


Research conducted during the past year under Contract No. N00014-78-C-0889 
may be conveniently divided into four categories. 

I. Baroclinic Tsunami Generation . Utilizing experimental data obtained 
previously, a study was completed on the internal waves generated in a stratified 
ocean by a vertical motion of the sea floor. Details of this study are presented 
in the enclosed manuscript (see Appendix II) which will appear in the Journal of 
Physical Oceanography. 

II. Evolution of Long Internal Waves . Subsequent to the generation processes 
mentioned above, the evolution of long internal waves on a thin pycnocline has 

been examined. Dr. Harvey Segur of Aeronautical Research Associates of Princeton 
and I have analyzed the experimental data and prepared a draft manuscript which 
examines several model equations for internal wave evolution in the context of 
experimental data. This work is ongoing; a copy of our final manuscript will 
be forwarded as soon as possible. 

III. Long Waves Generated by a Complex Bottom Motion . As a final statement 
on barotropic water-wave generation by motions of the sea floor, Professor Frederic 
Raichlen (of Cal Tech) and I completed a manuscript (enclosed) which illustrates 

a technique for computing wave structure when the time-displacement history of 
the sea floor is very complicated. The analytical mode was verified by experiments. 
These results were presented at the 17th International Conference on Coastal 
Engineering in Sydney, Australia (March 1980), and the manuscript (see Appendix III) 
will appear in the meeting proceedings. 






IV. Lecture Notes on Small-Scale Ocean Waves. At the invitation 


of the Italian Society of Physcis, I presented a series of lectures at their 
summer course on "Topics in Ocean Physics" held during 7 July - 19 July 1980 at 
the International School of Physics "Enrico Fermi", in Varenna, Italy. These 
lectures focused on the nonlinear aspects of water waves including solitons and 
resonant-interaction forcing of edge waves. The collection of lecture notes 
from the ten invited speakers will be published by the North Holland Publishing 
Company. A copy of my notes is presented in Appendix IV. 

V. Evolution of Short-Wave (Envelope) Solitons . The first goal of the 
experimental research program on short-wave evolution was the direct generation 
of an envelope soliton. The best results obtained to date for a wave packet 
with a carrier wave period T = 1 sec and water depth h = 1 m are shown in figures 
1 abed. Results at four positions along the (Berkeley) wave tank are presented 
with theoretical soliton profiles based on the local amplitude of the wave packet. 
The initial data measured at x = 6 m from the wave paddfe agree well with the 
theoretical soliton profile. During subsequent propagation the packer retains 
its compactness and travels with the (predicted) linear group speed of the carrier 
wave. However, the measured packet deviates slightly from the predicted soliton 
profile; this behavior is especially obvious at the last station (x = 45 m) of 
measurement. The front-to-back asymmetry emerging at x = 45 m is similar to that 
observed by previous researchers on deep water wave packets. In fact, the 
symmetrical measurements prior to the last measurement station are quite exceptional 
and have not been reported previously in the literature. (I might add that the 
first three wave traces of figure 1 are almost identical to those generated in 
the Florida facility; only the added length of the Berkeley wave tank permitted 
the distortion to be clearly observed.) 
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Below each wave trace in figure 1 I have also presented a histogram 
of the node-to-node wave periods through the wave packet. This rather 
crude attempt at complex demodulation provides some insight into the wave 
packet evolution. At the first two stations, local wave periods differ 
little from the programmed period of 1 sec. However, at the last two 
stations the leading waves are clearly developing periods up to 15% longer 
than the 1 sec period of waves near the center of the packet. The origin 
of the longer wave content is not clear although its emergence at the 
front of the packet is probably responsible for the changes between meas¬ 
ured and predicted packet piofiles. The most plausible explanation for 
the unwanted wave content of the packet is that it is a consequence of 
imperfect generation by the wave maker. It is an unfortunate (theoretical 
and observed) property of short wave solitons that these nonlinear waves 
propagate with the linear group speeds of their carrier waves. Hence, 
unlike their long-wave counterparts, short-wave solitons cannot separate 
from any background "noise" with the same frequency. (The background 
noise corresponds to the radiation or non-soliton component of the solution 
to the nonlinear Schrodinger equation.) The mergence and slow migration 
(at linear speeds) of the longer wave components in figure 1 may he a 
i ,uiifestation of this phenomenon. In this regard it is important to note 
that the radiation wave components remain dispersive even at nonlinear 
order and their amplitudes should decay in time t with an inviscid decay 

-h 

rate of t . Since the inviscid soliton does not decay, when noise is 
present one expects the measured data to agree hotter with the theoretical 
profile with distance (or time) down the tank. fxactly the opposite 
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bebavi or is observed in figure 1; however, it mist be remembered that glow¬ 
ing distortions might still occur until the radiation attains its asymp¬ 
totic character. 

Although inviseid solitons do not decay with time, water-wave solitens 
are subject to viscous attenuation of their amplitude. Two models for the 
viscous damping of short-wave solitons have been developed independently by 
Dr. Harvey Segur and myself (with considerable assistance from N. K. Lin, 
a graduate student at Berkeley). Both models predict similar results even 
though the approaches to model development were quite different. Experi¬ 
mental verification of the models must be considered limited for large pro¬ 
pagation distances due to the distortion discussed previously; short distance 
evolution according to these models is verified by the measured data. As 
in all \i •" mis damping models, contamination of the water surface enters as 
an e pi l ie. a 1 loefficient ranging from zero for no contamination to unity 
for a fully contaminated (and immobile) surface. Theoretical results for 
these two limiting cases are shown in figure 2 with experimental data for 
the wave evolution in figure 1. Measurements suggest that some surface 
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c on t a i n a t i on is present in the Berkeley data which is expected since water 
:iins in the tank continuously. Theoretical and experimental results 
for the ‘-.esc initial data observed in the Florida facility are also shown 
in figure 2. In Florida, well water was pumped to the wave tank immediately 
prior to the experiments; lance, surface contamination should he less; this 


expectation is verified by the data. (Ihc narrower tank in Florida leads 
to greater wave attenuation with distance of propagation.) 

Tt should also be noted that the viscous damping models for envelope 
solitons indicate that they decay twice as fast as a linear^ wave train of 
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the same frequency. Unfortunately, water wave trains of sufficiently small 
amplitude were not measured to test this prediction. However, measurements 
of wave trains with the same degree of nonlinearity as the comparison soli- 
ton (prior to the visual onset of Benjamin-Feir instability) indicate that 
the solitons decay only one-half as rapidly as the uniform wave train. 

The evolution of arbitrary packets of initial data was also investi¬ 
gated. An example of results is shown in figures 3a-b at two measurement 
stations; histograms of wave periods are shown in figures 4a-b for each 
measurement. The dominant wave packet in the downstream measurement of 
figure 3b clearly has evolved a soliton profile. The remaining wave packets 
do not have soliton profiles and represent the radiation component in the 
initial data. The period histogram of figure 4b clearly shows the linear 
dispersive nature of the constituent wave components which are rank ordered 
by their group speed. (Also see Appendix IV). 

Experiments on the stability of 2-D envelope solitons to 3-D perturba¬ 
tions have proven to be inconclusive. Very small 3-D effects were measured 
in some experiments but without a consistency or magnitude sufficient to 
infer definitive behavior patterns. The experiments did highlight a defi¬ 
ciency in the analytical studies on this aspect of envelope soliton. All 
explicit criteria which provide definitive predictions on the tank (or 
crest) widths necessary for instability assume infinite water depths. How¬ 
ever, the nonlinear description of these waves contains terms which decay 
algebraically--not exponentially--with depth. Hence, it is dangerous to 
apply these criteria to finite depth experiments, even for kh ~ 10. 
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Finally, based on the data presented here as well as many other experi¬ 
ments on varying aspects of short-wave solitons, it is clear to me that the 
reality of these systems is more difficult to disentangle than that of their 
long-wave counterparts. However, even with the confusion arising from 
microscopic observation of these waves, it is also apparent that the macro¬ 
scopic structure predicted by the nonlinear Schrodinger equation is present. 
All of these results emphasize the need for more analytical study coupled 
with delicate experiments capable of yielding quantitative measurement of 
third order (and higher) parameters. 
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Baroclinic Tsunami Generation 

Joseph L. Hammack 

Department of Civil Engineering, University of California. Berkeley 94720 
(Manuscript received 13 August 1979. in final form 2 May 19K0i 

ABSTRACT 

An analytical and experimental study of the baroclinic waves generated bt a monopole dislocation 
of the sea floor is presented. Analytical results are based on a two-dimensional and linearized description 
of motion using a two-layer approximation for density variation; experiments utilize a stratification 
with finite (nonzero) pycnocline thickness. Scaling parameters which characterize the generation process 
and the potential role of nonlinear effects are discussed. It is shown that the barotropic modes are 
not affected by the small differences in potential density typical of ocean stratifications and all previous 
results for these waves are applicable. The two-layer approximation is found to provide an accurate repre¬ 
sentation of the (long) baroclinic waves typical of tsunamis. Like the barotropic response, baroclinic 
generation is impulsive and linear resulting in wave amplitudes proportional to the vertical offset of the 
sea floor. Near the generation region barotropic waves have amplitudes of one-half the sea floor displace¬ 
ment while the baroclinic waves are attenuated further by the ratio hjh, where h is the total fluid depth and 
h, the upper layer thickness. Although Coriolis effects are not included in either the analytical or experi¬ 
mental models, these effects may often be significant for baroclinic waves. In general, the potential role 
of Coriolis forces is both earthquake and site specific Regardless, the analysis herein remains applicable 
for times smaller than the local inertial period 


1. Introduction and major conclusions 

Tsunami research to date has generally focused on 
the waves at the ocean's surface, ignoring vertical 
variations in potential density which typically exist. 
Submarine earthquakes which produce significant 
surface (barotropic) waves will also generate in¬ 
ternal (baroclinic) tsunamis in stratified regions. 
These baroclinic modes were examined first by 
Cherkesov (1968) who presented a theoretical in¬ 
vestigation of waves generated by a motion of the sea 
floor in an ocean with an exponential increase in po¬ 
tential density with depth. Subsequently. Fedo- 
senko and Cherkesov (1968) presented similar re¬ 
sults for a two-laver stratification which is more 
representative of actual ocean conditions. Solutions 
were based on a linearized description of motion for 
both two-dimensional and three-dimensional (axi- 
symmetric) source regions; Coriolis forces were 
neglected. In both presentations, the integral solu¬ 
tions were evaluated only for the far-field at large 
distances from the source region. These results 
showed that a displacement of the sea floor pro¬ 
duced (fast) barotropjc waves and (slow) baro¬ 
clinic waves which separated with time. Each wave 
system induced a coupled motion at the opposite 
interface. It was further claimed that the asymp¬ 
totic amplitudes of the internal waves could exceed 
those at the free surface by a factor of 2? for 
ocean conditions. Neither the relevance of asymp¬ 
totic linear analysis to actual tsunamis nor any 
limitation of the results due to omission of Coriolis 
forces were discussed by the authors. [Hammack 
and Segur (1978) have shown that linear asymp¬ 
totic solutions never become applicable for baro¬ 
tropic tsunamis.) 
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The objective of the present paper is to examine 
the near-field behavior of baroclinic waves gener¬ 
ated by vertical motions of the sea floor. Both 
analytical and experimental results are presented 
for a simple model of generation consisting of a 
block upthrust (or-downthrow) according to a pre¬ 
scribed time-displacement history. Analytical re¬ 
sults are based on a two-layer stratification of 
potential density while the experiments utilize a 
stratification with finite (nonzero) pycnocline 
thickness. 

The organization of this paper and some of the 
major results are as follows. In Section 2 a brief 
description of tsunamigenic earthquakes is pre¬ 
sented and the model to be examined herein is 
introduced. Nondimensional parameters which 
characterize the generation process are discussed 
in Section 3 where it is shown that the generation 
process of actual tsunamis is impulsive and linear 
for both the barolropic and baroclinic response. In 
Section 4 the linear solution for the baroclinic 
model is developed and other well-known properties 
of long internal waves at the interlace of two 
liquids are reviewed. In the Boussinesq limit of 
small density differences, it is shown that the baro- 
tropic response is not affected by the stratification. 
In addition, both barotropic and baroclinic tsunamis 
have initial shapes similar to the permanent de¬ 
formation of the sea floor. The maximum amplitude 
of the surface wave is one-half the vertical dis¬ 
placement of the sea floor while the maximum 
amplitude of the internal wave is additionally 
attenuated by the ratio of upper layer depth to total 
depth Experimental equipment and procedures are 
described in Section 5 and theoretical results are 
compared with the experimental measurements in 
Section 6. Results indicate that the maximum 
amplitudes of the barotropic and baroclinic waves 
near a source region may be represented by a single 
(and simple) functional relation over the full range of 
their respective generation time scales (see Fig. 10). 
The potential role of Coriolis effects for baroclinic 
waves is discussed in Section 7, where it is dem- 
onstiated that time scales of baroclinic modes can 
a 1 nui) often approach the inertial period. 


2. Tsunamigenic earthquakes and the proposed 
model 

Most large tsunamis are associated with sub- 
oceanic earthquakes, centered near plate margins, 
during which significant vertical offset of the sea 
floor occurs. In order to determine the waves 
generated by mechanisms of this type, a knowledge 
of the water -depths, stratification and latitude at 
the source site as well as the spatial and tem¬ 
poral distributions of the sea floor dislocation is 
required. Since the objective herein is to delineate 
fundamental relations between the generation 
mechanism and wave structure, we seek a tractable 
model which embodies essential source parameters 
in a simple way. With this goal in mind we will 
assume a uniform depth It in the generation re¬ 
gion and adopt a two-layer model to represent 
the ambient stratification in potential density. In 
addition we will neglect Coriolis forces (and their 
latitude dependence); limitations of the results due 
to this omission will be discussed a posteriori. 
Before introducing further simplifications, a brief 
description of features common to tsunamigenic 
earthquakes and their (dimensional) scales is 
presented. 

It is both convenient and appropriate to define the 
source size of the tsunami (and earthquake) as 
the area affected by aftershocks immediately 
following the generative earthquake. This area is 
typically elongated and elliptical in shape with the 
ruptured portion of the fault joining the foci of the 
ellipse. Early efforts to correlate the source area 
S with Richter magnitude 1 M s proved unsatisfactory 
due to the marked scatter in the data. Kanamori 
and Anderson (1975) inter alios have shown that 
definitive correlations exist between S (km 2 ) and 
the seismic moment M n (dyn cm). Representative 
values based on their results for interplate earth¬ 
quakes are shown in Table 1. In order to provide 
perspective, it is still desirable to relate the seismic 
moment to Richter magnitude. As expected, results 
here are less definitive and no single empirical 
relation exists; however, Kanamori and Anderson 
(1975) do present correlations for three classes of 
earthquake magnitudes. Approximate values of M, 
based on their results are shown in Table I. As a 
first approximation, we may assume that the source 
areas of Table 1 represent regions of monopole 
vertical offset (upthrust of downthrow). With this 
simplification the source size can be characterized 
by a single length dimension L ~ 5 1 ' 2 ; representative 
values are shown in Table I. If we further assume 
the static and seismic moments are equal, mag¬ 
nitudes of the average (and permanent) vertical 
offset £„ can be found by 
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!C,,| = A iJyS, 

where n is the shear modulus (rigidity) of the frac¬ 
tured medium. Estimates ot |£„| with p. = 5 x 10" 
dynes cm 2 are presented in Table 1. The temporal 
history of the sea floor dislocation is characterized 
by a minimum of two time scales: a time t c cor¬ 
responding to the time required for the vertical 
offset to occur at a specific spatial location and a 
time t„ representing the period required for the 
rupture to propagate along the fault. Empirical 
results of Kanamori and Anderson (1975) suggests 
that t H ~ 10/ r for tsunamigenic earthquakes (M s 
s 6.5) with estimates of the rupture time ranging 
from t R - 10 s for moderate earthquakes to t R 
~ 100 s for earthquakes with fault lengths exceeding 
100 km (based on rupture speeds V R ~ Llt R ~ 2-3 
km s' 1 ). As suggested previously, vertical offset is 
primarily responsible for tsunami generation; hence, 
the time t r is most relevant for subsequent 
scaling. 

In light of the above description, we will adopt 
the following simplifications for the generation 
model. First, the smallness of the linear source 
dimensions of Table 1 relative to the earth s 
radius permits a plane-earth approximation. Second, 
for monopole and uniform vertical offset, the 
elongation of the source dimension suggests that 
initial propagation of waves near the center of the 
source region is one-dimensional and normal to the 
fault. Hence, we arrive at the simple generation 
model shown in^ Fig. I. The actual sea floor dis¬ 
placement used in computations is given by 

£(.r./) = ^„(e - 1 )H{b 1 - * 2 ), (1) 

where H represents the Heavyside stepfunction. In 
the dislocation of (I), we have taken 2b ~ L and, for 
experimental convenience, we will take 



a 


3. Scaling the generation process and a comment on 
nonlinearity 

In order to provide a rational basis for choosing 
computational examples, nondimensional param¬ 
eters which characterize generation and typical 
magnitudes for tsunamigenic earthquakes are re¬ 
quired. With this objective in mind, a brief review 
of results by Hammack (1973) for generation of 
barotropic waves is useful. Three nondimensional 
parameters are required to characterize barotropic 
generation for simple models similar to that dis¬ 
cussed in Section 2: an amplitude scale. |(„| Ih: a 
size scale, bllr. and r = t r C,/b termed a time- 
size ratio where C, = (gh ) 12 is the speed of long 
barotropic modes. [Alternatively, r may be con¬ 
sidered an inverse Froude number analogous to that 
used by Noda (1970) for landslide-generated waves.) 
The generation process may be classified according 
to the time-size ratio as impulsive for t 4 I, creep¬ 
ing for t > 1 and transitional for r = 1 . The specific 
range of t for each classification is only a weak 
function of the size scale with the transitional regime 
essentially disappearing for large size scales (b/h 
s> 1). The amplitudes of waves propagating out of 
the generation region for an impulsive motion are 
proportional to the bed displacement (rj cc („) 
where the constant of proportionality depends 
only on the size scale. For b/h > I, this constant 
attains its maximum value of one-half and the 
tsunami near the source region is similar in shape to 
the permanent deformation of the sea floor. For 
creeping bed motions, wave amplitudes are reduced 
in inverse proportion to the magnitude of the time- 
size ratio (17 °c £ 0 It). Nonlinear effects become 
significant when generation is impulsive and 
3 0.2. Assuming depths in the generation region 
of h = 200 m on the shelf and h *= 4000 m in deep 
water, size scales for tsunamigenic earthquakes are 
typically large (based on the lengths L ~ b given 
in Table 1); hence, energy is concentrated at long 
wavelengths. In addition, amplitude scales are small 
based on values of |£„| given in Table I. Finally, 
using characteristic times up to t,. — 10 s. the baro¬ 
tropic response is typically impulsive (- <S 1). Based 
on these characteristic values, the barotropic re¬ 
sponse is expected to be linear and the wave struc¬ 
ture near the source region should resemble the 
permanent deformation of the sea floor with a 
50% attenuation in amplitude. It is noteworthy, and 
fortunate, that the barotropic response is repre¬ 
sented by the simplest possible limiting conditions 
of the generation parameters. 


Physically. t c represents the time required for two- 
thuds of the vertical offset to occur. 
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An examination of the independent variables 
(p 1 ,p 2 ,/i,,A 2 ,A,r,.,£ u ,g) for generation in a two-layer 
ocean (see Fig. I) indicates that five nondimen- 
sional parameters are required to characterize 
generation. A possible parameter choice for small 
density differences is 

blh 

7, = t,.CJb = Inf’ll)'Vb , (3) 

r-j = i r C t lh = IAeghJiJh)'-/b I 

(/ t ,// I ,)' 1 j 

where h - /i, + /i-> and e = (p 2 - Pi)/p-> 1. Of 

these five parameters, the first four preserve the 
scaling characteristics of the barotropic mode and 
indicate the expectation of a similar role for a time- 
size ratio (r-i) based on the internal (long) wave 
speed C> (see Section 4). Notice that the two 
time-size ratios satisfy t-Jt , = 0 (e' 2 ); hence, bed 
motions which are classified creeping for the surface 
response (t, ^ 1 ) may be impulsive for the inter¬ 
facial response (r 2 1). Since r, < 1 for prototype 

tsunamis, the baroclinic response is always impulsive. 

The fifth parameter |(/i,//i 2 ) - 1) is geometric in 
nature; the various roles of this parameter and its 
related forms A,//r or /;.,//( will appear in Sections 
4 and 6 . Its specific role as an indicator of the 
strength of nonlinear effects during generation can 
be established as follows. In a two-layer system, 
nonlinear effects may occur for both the surface 
and interfacial motion. For the free motion of long 
internal waves in a two-layer system with small 
density differences and uniform depth, it can be 
shown (see Keulegan, 1953) that the appropriate 
measure of small but finite nonlinearity is |t) 2 (/i, 
- A 2 )//i,A 2 |. It will be demonstrated in Section 4 
that t), cc £ 0 (/i,//i) for impulsive generation; hence, 
nonlinearity is measured by the parameter |(£ 0 / 
li)Ui t ili-> - 1 )|. Since hjh> =£ l for typical ocean 
stratifications and (£ 0 //i) is small, both barotropic 
and baroclinic free motions are linear initially. 

4. Theoretical analysis 

We consider a two-dimensional ocean as shown 
in Fig. I consisting of two fluid domains D x and 
D, bounded above by a free surface S,, below by a 
solid boundary 5* and contacting at a liquid-liquid 
interface 5,. Both fluid regions are unbounded in 
the direction a of wave propagation. The upper and 
lower fluids have densities p, and p 2 , respectively, 
and are gravitationally stable (p, < p 2 ). Initially, the 
two fluids are at rest with uniform depths li, and 
A,. Subsequently (/ >0), the solid boundary S,, 
is permitted to move in a ptescribed manner given 
by y = -bn + ((a,/), where £(*,/) is given by ( 1 ). 
The resulting motion of the free surface ry,(.v./) and 
interface T? 2 (.r,t) are required. Since both motions 
are probably linear for prototype tsunamis, a first- 
order (linear) description of motion is employed 
from the outset. 
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Assuming the fluids in each domain to be incom¬ 
pressible and their motion irrotational, velocity- 
potentials 4>,(x.y j) and <j>.,(x,y,t) exist for Z), and 
D... respectively. Conservation of mass requires 

r-<i>, =0 in D,, ( 4 a) 

V= 0 in ( 4 b) 

where V- is the Laplacian operator. Assuming the 
fluids to be inviscid, the kinematic boundary condi¬ 
tions at each surface are 

</>,. = T) ti on y = A,. (5a) 

= V>, on y = 0. (5b) 

By further neglecting surface energy at S f , the condi¬ 
tion for constancy of pressure across and along the 
free surface is 

4>i, + 8Vt = ° on y = /i,. ( 6 a) 

(The hydrostatic pressure p,g/t, has been adsorbed 
into <bi.) Continuity of pressure across 5 , requires 

+ £92 = — (<£,, + gr) 2 ) on y = 0 . ( 6 b) 

P 2 

The kinematic condition at the solid boundary is 

4>i, = £< on y = —/i 2 , (7) 

where it is assumed that motion of the solid boundary 
is not affected by fluid motion. 

Combining the kinematic and dynamic boundary 
conditions at S f and S t yields 

<*i„ + 8<t>i, = 0 on y = A,, ( 8 ) 

^ 2 ,, + = — (d>,„ + g<t>, ) on v = 0. (9) 

Pi 

Using the Laplace transform in t and the Fourier 
transform in v. Eqs. (4), ( 8 ), (9) and (7) become 


^i„ _ k-4>, = 0 in £>,, (10a) 

- k'4>i = o in Z? 2 , (10b) 

s ‘<t >i + = 0 on v = A,, (10c) 

.t 2 <f>, + £</>,„ = — (s'-'<j>, + gi>, ) on v = 0, (lOd) 

Pi 

4>i, = si on y = -h 2 , (lOei 


where the tilde of a function/(v.r) indica'e* 
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/(A,.v) = I dx ( e >kJ ’e~’f(x,t)dt. (11) 

The function/U./) is recovered from (11) by invert¬ 
ing the Laplace and Fourier transforms, i.e.. 

1 P n 

./(*./> = — dk 
27T L* 

lim — f ' *- kT " s 'j\k,s)ds , (12) 

' 2m 

using the complex inversion integral for the Laplace 
transform. The transformed surface and interfacial 
waves are 



V\ = 


on v = //,. 


(13a) 


Vi = — l(Pi/p*)«^i - </>*! on v = 0. (13b) 

«.t? 

Solution of (10a,b) yields 

<i>,(A.v,.v) = A, coshAy + A 2 sinhAv, (14a) 

<h,(k,y,s) = A-, coshAy 4- A, sinhAv, (14b) 

and the coefficients A, = A,(A.s) with i = 1,4 are 
found from Eqs. OOc)-(lOe) to be 


A i = 


" O 


A, - A, 
A i 


Dk coshA/i , 

-v :, A 


•v-7-, + ,</A 1, 


(IA coshA/n, 


I v s + cA7, 


(/A coshA/i, 


.v‘p r 7, + s-gk + e( jt?A _ ) 2 7~,). 


(15a) 


(15b) 

(15c) 


with 


8 = s 4 ( I + — 7,7,) + v 2 (7, + 7,)i>A 

’ />2 ' 

+ t(gA)-r,T,, (16a) 

T, = tanhAA,, i = 1.2. (16b) 

Finally, the transformed free surface and internal 
waves are found from (13) to be 


a 4 

(17a) 

tl coshA/i i coshA/z 2 

- - IP oWI. 

0 coshkh-z 

(17b) 


The transform of the bed motion given by (1) is 


aA-.v) = 2£ 0 


sinAA> 


a 


k Li(s + a) 


( 18 ) 


Substituting (18) into (17), performing the integra¬ 
tion around the Bromwich contour, taking only the 
real part and noting that the integrands are even 
functions of A, the free surface and interfacial 
waves are found to be 


T) ,(*./) 


T /•(*./) 


r B,Bi 
coshA/i, coshA/t., 

x [e-“'- B.,B, + B :r B 6 ]dk, (19a) 


B.B, 


coshA/to 


i + 



5.5 4 


with 


B i 

« 2 

S, 

5.V 




</A. 


(19b) 


-2/, sinA/i cosA.v 
;rA[l + <p,/p 2 )7,7 2 ] 


(20a) 


«■* 

(or + ci> t 2 )(or + (a, 3 ) 
ai, 3 ( or + ta 2 2 ) 
of(ai| 2 - o>> 3 ) 

1D| . 

COSoj,/ + — Sinca,/, 
« 

ta/fu 3 + a>, 2 ) 
or 2 (<o,-’ + <a 2 2 ) 

cosw 2 / +■ — sima./, 
a 


(20b) 

(20c) 

(20d) 

(20e) 

(200 


and u>,- U = I, 2) are the roots of the dispersion 
relation 

(1 + p,/p,7 ,T 2 )oj 4 - gk(T, + T 2 )oj- 


given by 


+ e{gktt,T, =0. (21) 


! (7-, + T,) ± |(T, + 7 2 ) 2 

, j - 4e(l + p,/p 2 7,7 2 )7,7 2 ]'' 2 

<uT. 2 = gA j - . (22) 

2(1 + p,/p 2 7,7 2 ) 

(Note that io = -is in the above analysis.) Eq. (21) is 
well known for two-layer systems with a,, 2 corre 
sponding to the barotropic modes and o> 2 2 corre 
sponding to the baroclinic modes (see e g I n"4 
1932 p 37?) 
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Although the Fourier integral solutions given by 
(19) are algebraically complex and the,wave struc¬ 
ture is difficult to see, it is easily shown that these 
solutions contain more familiar results. For ex¬ 
ample, in the Boussinesq limit of small density 


differences <e 

—* 0) we find 




1 ‘ cosA.v sinAh 

f a " ) 


Vi V - 

7 r . 

A coshA/i 

1 <r + oj- ' 


X 

( (- - COSOlf 

- — sincof )<(A. 

(23) 


to," ~ or == f>k 

tanhA/i, 

(24) 


where 17 is the solution found by Hammack (1973) 
for a barotropic ocean. Hence, the surface waves 
in a two-layer system are not affected by the strati¬ 
fication in the Boussinesq limit. Since density dif¬ 
ferences are typically small for oceans U ~ 0 . 002 ), 
previous results for surface waves in a barotropic 
ocean are applicable to leading order. 

It is also useful to approximate specific properties 
of the internal and surface waves for long modes 
such that 


e a = O(A/t,) = O(A/ift <§ 1. (25) 

P-.- 

Under these conditions the well-known phase 
speeds of the surface and internal waves are found 
(from Eq. ( 22 )] to be, respectively. 


= 1 or Ik- - git 1 


It 


— + CKe 3 ) 
It, 


(26a) 


CV = ~ 



0 (£-')|. 


(26b) 


The significant difference in phase speeds for these 
long waves. CJC, = Ole 1 ' 4 ), suggests that even 
though the surface and internal waves overlap 
spatially near the generation region, they will 
separate 1 ipidly and nonlinear interactions between 
(he two systems will not have time to become sig¬ 
nificant 

In addition to the phase speeds it is also of interest 
to examine the amplitude of the motion induced by 
each wave system at the opposite interface. Follow¬ 
ing Keulegan (1953), normal-mode analysis sug¬ 
gests that 

ft., = — .V,, (27a) 

It 


ft, = 


€ 



where 0., is the coupled motion of the interface due 
to a wave /V, at the free surface and ft, is the coupled 
motion of the free surface due to an internal wave 
.Vj. Here, the waves ,V, and /V, represent the inde¬ 
pendent modes at the free surface and interface, 
respectively. (Note that ,V, = 17 , and N 2 = r), when 
the surface and internal waves do not overlap 
spatially.) The motion ft 2 is barotropic and would 
exist at the interface location even if the fluid was 
not stratified. 

In order to see the wave structure in (19), numeri¬ 
cal solutions at r = b for both barotropic and baro- 
clinic waves propagating from the source region are 
shown in Fig. 2. Results are presented for impul¬ 
sive response at the free surface and interface 
(Fig. 2 a) and for impulsive interfacial response with 
creeping free surface response (Fig. 2b). Magnitudes 
of the nondimensional generation parameters used in 
computations for Fig. 2a are chosen to be asymp¬ 
totically representative of prototype phenomena: 
dimensional quantities correspond to experiments 
which will be discussed in Section 5. 

As expected, the surface waves do not appear to 
vary with layer-depth ratio in Fig. 2. In fact these 
results are identical to those for a homogeneous 
ocean computed by (23) except for a small residual 
lowering of the water level which persis*s until the 
internal wave passes. It should be emphasized that 
surface and internal waves overlap spatially in the 
generation region; hence, these results represent a 
linear combination of both the independent (,V) and 
induced motions (0) at the interface and free 
surface. The residual lowering of the free surface 
is the 0(e) motion /3, [see Eq. (27b)] induced 
by the internal wave. In fact, this induced motion 
differs for each depth ratio; however, the small 
differences are not discernable at figure scale. The 
impact of the surface wave at the interface, 
ft ; , is more pronounced and clearly evident in Fig. 
2a. This induced motion is easily estimated from 
(27a) using the free surface solution [since iV, = 77 , 

+ O(e)]. When ft 2 is eliminated from the internal 
wave computations, iV 2 results and is also shown in 
Fig. 2a. The leading waves of the independent 
internal and surface modes are similar to the per¬ 
manent deformation of the sea floor. The maximum 
amplitude of the surface wave is one-half the sea 
floor uplift, i.e., (,V,) inax = £„/ 2 , as found previously 
for a homogeneous ocean. The maximum amplitude 
of the independent internal waves (/V.,)„,.„ for these 
impulsive motions appears to vary with the layer- 
depth ratio according to 

GVftmax = <Mt)GV P >max = (£„/2)( It ,//l). (28) 


(27b) 






JPO—37201 Sept' 


Gal 148 


This important result and further clarification of the 
effects produced by varying the size scale blh as well 
as relative thickness of the upper layer may be estab¬ 
lished in the following manner. Restricting atten¬ 
tion to (mathematically) impulsive bed motions, 
the initial (t = 0 + ) amplitudes at x = b can be eval¬ 
uated in closed form to yield in the Boussinesq 
limit 


2 £„ , rr b i 

r) i{h.O') = - arctan jtanh — — J , 


(29) 


T?,4/>.0*) 


Co 

= — arctan 


sinh( 77 />//t)/cos| ~ ) 


(30) 


The result for the surface wave (29) is identical to 
that for a homogeneous ocean (Sells, 1965) and 
clearly shows that tj, — N, — 17 — Co/2 for blh > 1. 
Eqs. (29) and (30) are shown in Fig. 3 as a function 
of blh and hjh. (Note that (29) and (30) are equiva¬ 
lent when h,lh —> 0.] Results for the internal wave 
are identical to the surface wave for blh > 1 
regardless of the layer-depth ratio. For/;//; < 1, the 
internal wave amplitude exceeds that at the free 
surface as hjh increases (e.g., hjh > 0.25). In nil 
cases, the amplitudes become inversely proportional 
to the size scale as blh decreases. Noting that 
iV -,(b.0 *)' — T] 2 (b,0 r ) - f} 2 {b,0 + ), where p is given 
by (27a), Eqs. (29) and (30) may be used to deter¬ 
mine the behavior of the internal mode with size 
scale; the result is 


O') 

= k i 

• n j 

i 


arctan 


sinh(7r/;//i) 


cos 


| 

2 h ) 



x arctan 



(31) 


For blh > 1, Eq. (31) yields, in closed form, the 
numerical results found previously for the internal 
wave amplitude given by (28). Hence, even though 
the free surface and interface are lifted (lowered) 
the same amount by an impulsive bed motion with 
blh > I. part of the potential energy gained at the 
interface is associated with the barotropic motion 
P<. Consequently, the amplitude of the internal 
wave is less (by the factor hjh) than that of the 
surface wave where all of the potential energy 
gained, to O(e), is available to the barotropic mode. 


5. Experimental equipment and procedure 

A seriesofexperiments was conducted in a labora¬ 
tory wave tank which is approximately 32 m long, 
60 cm deep and 40 cm wide The tank is equipped 
at one extremity with a wave generator consisting of 
a movable bed section (piston) whose time-dis¬ 
placement history is controlled by an electro- 
hydraulic-scrvo system. The wave generator can 
accurately model the bed motion described by ( 1 ). 
Both the wave generator and tank have been 
described in detail by Hammack (1972). 
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In order to model the density stratification of Fig. 
1, freshwater and brine are used. The miscibility 
of these fluids produces a finite pycnocline thick¬ 
ness 6 ; thus, this stratification provides a better test 
of the two-layer model as a useful approximation for 
thin pycnocline regions. The tank is stratified by 
first introducing a layer of freshwater with a depth 
/' 1 - Brine is then introduced beneath the fresh¬ 
water at ports distributed along the tank bed until 
a total depth h is attained. The volume of the 
reservoir used to mix the brine limits the depth It, 
to a maximum of 10 cm. For the experiments 
presented herein, the total depth h is fixed at 15 cm 
while the depth of the brine is varied with h< = 5 . 
7.5 and 10 cm. A density difference of ~5% is 
utilized in all experiments. This difference is suf¬ 
ficiently large to provide stability of the stratifica¬ 
tion during filling yet sufficientlysmall to provide an 
adequate test of analytical results based on small 
density differences. The length b of the movable bed 
section is constant (61 cm) for all experiments. 
Hence, the experimental size scale, blh = 4.07, is 
sufficiently large for results to be typical of the 
tsunamigenic case. Density stratifications are meas¬ 
ured using a conductivity probe at two locations 
before and after each series of experiments; typi¬ 
cal results are presented in Fig. 4 where /i, = 5 cm 
and li 2 = 10 cm. The similarity of initial stratifica¬ 
tions at both locations demonstrates the uniformity 
of this stratification along the wave tank. A thick¬ 
ness 8 for the pycnocline region is calculated 
based on the definition 



yielding initial values of 8 0 = 0.9 and 1.0 cm and 
final values of 8 , = 2.1 and 1.8 cm for the upstream 
and downstream positions, respectively. Actual 
stratifications are qualitatively similar to theoretical 
solutions of the diffusion equation for an initial 
density profile with a simple discontinuity at one 
level. Using this solution {p = p, 4 - dp erfc[y/ 
(4Dt) , J J} with the initial discontinuity occurring 
at y = 0 , the characteristic thickness 8 defined by 
(32) is related to the diffusion coefficient D by 


8 = (4-nDt)' -. ( 33 ) 

Hence, an apparent diffusion coefficient can be cal¬ 
culated for the time interval of the experiments by 


I / (S,) 2 ~ (8„) ; 

4n \ At 


( 34 ) 
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The apparent diffusion coefficients for the measure¬ 
ments of Fig. 4 are D = 2.3 x IO~-> cm 2 s'* at the 
upstream position and D = 1.7 x 10 -5 cm 2 s 1 at the 
downstream position. The more intense mixing at 
x = b results from the shedding of a vortex at the 
edge of the bed section during its motion (see 
Hammack, 1972). The mixing of the interface at this 
position was visually detected by adding blue dye to 
the brine. At the downstream position the apparent 
diffusion coefficient is only slightly greater than 
the molecular diffusion coefficient (D,„ = 1.5 x 10 ~ 5 
cm 2 s' 1 ) for salt in water. Hence, the shear layer 
developed during the passage of internal waves is 
laminar and does not mix the fluids appreciably. 
Each series of experiments was stopped when the 
interfacial thickness at x = b reached —2 cm. 

In order to quantitatively measure the internal 
waves, a laser-optics detector system is used. A 
schematic drawing of the internal wave gage is 
shown in Fig. 5. The light beam from a helium- 
neon gas laser (wavelength = 632.8 nm with 0.5 mW 
output) is aligned perpendicular to the glass side- 
walls of the wave tank. A system of cylindrical 
lenses transforms the circular beam into a sheet of 
light with a constant vertical dimension of 4.5 cm 
which is directed horizontally across the wave tank. 
After traversing the tank, the light is focused onto 
a photodiode which provides an output voltage pro¬ 
portional to the incident light intensity. By dying the 
brine dark blue and permitting the interface to inter¬ 
cept the light sheet when the fluids are quiescent, 
subsequent motions of the interface induce changes 
in the output signal of the photodiode which are 
displayed on an oscillograph recorder. It should be 
noted that the dyed interface seen by the internal 
wave gage appears in the upper region of the dif¬ 
fuse salinity interface as indicated in Fig. 4. The 
laser, optics and detector system are mounted on a 
continuous frame shaped in an inverted U over the 
wave tank; hence, the wave gage is calibrated by 
raising and lowering the frame differing amounts 
during quiescent fluid condition. Calibration curves 
are inherently nonlinear for this measurement sys¬ 
tem. Internal waves with maximum heights < 4.5 
cm are easily measured; however, large-amplitude 
motions are generally avoided during wave measure¬ 
ments at t = b due to the vortex-induced mixing 
at the position. The vortex generates a high- 
frequency signal at the photodiode which is super¬ 
posed on the signal corresponding to the internal 
wave. In addition, the background signal for quies¬ 
cent conditions may change appreciably before and 
after an experiment. Surface waves are electroni¬ 
cally recorded using conventional parallel-wire 
resistance probes and an oscillograph recorder. 
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6 . Comparison of theoretical and experimental 
results 

In order to examine the near-field wave struc¬ 
ture, the temporal behavior at the edge of the gen¬ 
eration region (,r = b) is examined. Wave motions 
at this location are primarily free (unforced) and 
provide the initial condition for propagation models 
which may be patched to the generation model. 
Results are presented for a wide range of time-size 
ratios in order to more clearly delineate the tsunami- 
genic case (impulsive with It, « lu and bib > I). 

a. Wave profiles 

Comparisons of theoretical wave profiles at 
v = b with experimental data are shown in Figs. 6 
and 7. Results are shown for positive (uplift) and 
negative (downifyow) bed motions. Linear com¬ 
putations for the free surface motion agree well with 
measurements. (Surface wave measurements ate 
coincident for the two cases of bed uplift.) Linear 
results for the internal waves also agree well with 
measurements; however, there are differences 
which require explanation. Measured waves at the 
interface exhibit a high-frequency noise of small- 
amplitude superposed on the main wave. The source 
of this high-frequency noise is twofold. First, as 
mentioned in Section 5, a vortex is shed at the down¬ 
stream edge of the piston (.r = b) during rapid 
motions. The effect^of this vortex on the amplitude 
measurement at the interface is most pronounced for 
positive motions with thin lower layers (/i, < It,) 
and was visually observed to be responsible for 
the large depression in amplitude near Ug/h )' 2 = 12 
in measurement (d) of Fig. 6 . Second, gas bubbles 
formed in the thin gap between the piston and its 
confining tank walls. These bubbles were freed 
during bed motion and subsequently those from the 
gap along x = b penetrated the light sheet of the 
measurement transducer. (Brine was especially 
prone to bubble formation resulting from aeration 
during mixing and chemical reactions at unavoid 
able scratches in the anodized aluminum of the 
piston.) 
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The lime required for the lead wave to propagate 
past v = b differs from the linear prediction for 
both surface and internal waves in Fig. 6 . This 
behavior was observed previously for surface waves 
by Hammack (1973) and may be partially attributed 
to nonlinear effects arising from both finite wave 
amplitudes and finite changes in depth in the gen¬ 
eration region. Based on amplitude considerations 
alone, differences between linear analysis and 
measurements are expected to be maximum when 
/i, S h , and minimum when //, = h .. This expected 
behavior is exhibited in Fig. 6 even though ampli¬ 
tude nonlinearity is too small to account for all of the 
observed discrepancy. Another contributing factor 
is the eddy shed at x = h during piston motion. 
Fluid mixed by the eddy often remains in this re¬ 
gion for some time after the waves had departed. 
(Eventually, this fluid spreads into its layer of neutral 
buoyancy.) It was confirmed qualitatively by visual 
observation that this phenomenon causes an appar¬ 
ent increase in both the wave amplitude near its 
leeward portion and the lead wave period similar to 
the behavior observed in measurement (</) of Fig. 
6 (where /r, > h 2 ). 

Mode separation and motions induced at the op¬ 
posite interface are more clearly demonstrated by 
the downstream (.v = b + 20 h) results shown in 
Fig. 8 . Both experimental measurements and linear 
computations are presented for an impulsive bed 
motion with £ 0 = 1.5 cm,/i, = 5 cm and /i., = 10 cm. 
The surface wave arrives first at the downstream 
station at t(glh)' - = 20 corresponding to the speed 
C, given by (26a); the induced motion f3, at the 
surface is clearly evident. Observed and predicted 
behavior of the surface wave and its induced mo¬ 
tion agree well. Both the interface and free surface 
become quiescent prior to the arrival of the in¬ 
ternal wave. Predicted and measured shapes of the 
internal wave are similar; however, the maximum 
amplitude of the measured wave is 20 % less than 
predicted. In addition, the time required by the 
observed wave to reach the downstream station 
corresponds to an average speed which is 10 % less 
than the linear prediction for C 2 . These discrep¬ 
ancies probably result from the viscous and non¬ 
linear effects present in the experiments. 

h. Maximum amplitudes 

It is instructive to examine the variation in lead 
wave amplitude with the time-size ratio of the bed 
motion. In fact, results of this type provide the 
basis for classification of the generation process 
(impulsive, creeping and transitional). Fig. 9 shows 
experimental and theoretical results at v - b tor 
both the surface and internal waves. The maximum 
amplitudes of the independent modes (A/, and <V 2 ) 
are shown as a function of their respective time- 
size ratio and layer-depth ratio. Experimental meas¬ 
urements for both positive and negative bed motions 
are presented. (Nonlinear effects for all experiments 
shown in Fig. 9 are expected to be small based 
on the discussion in Section 3.) The time-size 
ratio for the surface response (r,) spans the 
classification range of impulsive to creeping while 
the internal response (t 2 ) ranges from impulsive 

I -'s'liu-' il 
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Surface wave amplitudes in Fig. 9 are identical for 
each layer-depth ratio (since the size scale bill is 
constant) and equivalent to that found for a homo¬ 
geneous ocean. Both experimental and analytical 
results in Fig. 9 demonstrate that CV,) nuis = £„/2 for 
impulsive motions while (.V, ) max a. £ 0 /r, in the 
creeping regime; previous results indicate that this 
behavior i? independent of size scale when bib > 1. 
The theoretical and experimental results for the 
internal wave amplitudes agree well. In addition, the 
variation of (/Vj) mitx /£n with r-< is similar to that of the 
surface waves except for an amplitude reduction by 
the factor hjh which appears to be valid even for 
time-size ratios exceeding the impulsive regime. 
This behavior permits the barotropic and baroclinic 
results to be collapsed into a single curve as shown 
in Fig. 10 for the full range of r, and r 2 (when the 
size scale satisfies bill > 1). The weak dependence 
(narrowing) of the transitional regime with larger 
size scales C blh 1) is also indicated in Fig. 10. 

7. A comment on Coriolis effects and the two-layer 
approximation 

It is customary and appropriate to neglect Coriolis 
effects for barotropic tsunamis since their energy is 
concentrated at periods (say. /, = 30-60 min) small 
compared to possible inertial periods. (Alterna¬ 
tively, the length scale L of the barotropic mode 
is small compared to the Kelvin-Rossby radius of 
deformation.) However, the sluggishness of the 
baroclinic response indicated by C 2 /C, — e 1 '* = 0.04 
immediately suggests the possibility of mesoscale 
periods for these modes. Since the length scale L 
for both barotropic and baroclinic modes is the 
same and fixed by the source mechanism, a first 
estimate of the dominant baroclinic period is 

h ~~ —— 1 , ~ «-"*/, ' 25/,. 

C 2 
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When/, is of order or greater than one inertial period 
U, = 12 h (sin<l>r' where <t> is the source latitude] 
the baroclinic modes must be considered as inertio- 
gravity waves. In fact, Coriolis effects will limit the 
period of these inertio-gravity waves to t ... < t, 
with t., —» i, as L —» (e.g., see Kraus, 1972, p. 

215). Neither the results of the preceding analysis 
nor the waveguide (one-dimensional propagation) 
approximation are applicable for describing the 
long-term (an inertial period on greater) dynamics 
ot these waves. However, the preceding results 
still remains applicable for propagation times much 
less than an inertial period. In either case, the 
variaodity of both source lengths for tsunamigenic 
earthquakes and the inertial period (with latitude) 
suggest that the potential role of Coriolis effects 
is earthquake and site specific. 

Finally, it should be noted that the two-layer 
stratification adopted herein permits an unbounded 
frequency spectrum for the baroclinic modes. In 
actuality, the continuous stratification of oceans 
will limit free baroclinic oscillations to frequencies 
less than the maximum Brunt-Vaisala frequency, 
say N,„. [As noted by Kraus (1972, p. 218), the two- 
layer model is equivalent to the first baroclinic 
mode of a continuous stratification with equivalent 
phase speeds.] For a typical ocean stratification 
with ,V,„ = 0.01 s ', baroclinic wave energy will be 
restricted to periods c, > 10 min. As already im¬ 
plied, energy generally will be concentrated at pe¬ 
riods much greater than this lower bound. 
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Fig. I. Fluid domain and sea floor dislocation model: 

(a) spatial, (b) temporal. 

Flu. 2. Theoretical wave profiles at x = b = 61 cm for« = 0.05 and h = 15 cm: 

(- 1 free surface waves (all It,lit): broken curves represent internal waves: 

l-) = 2/3: (-) It,lit = 1/2: (-) Mi = 1/3: I.) V.. (a) 

t , = 0.1. t. = 0.01: (b) r, = 5.2, t, = 0.5. 

Fig. 3. Theoretical Variation of maximum wave amplitudes with source size 
scale hilt at v = h for impulsive response: (solid line), (tj, ) ma ,/{o'. dashed line. 
( 1 1 ■: 

Fig. 4. Experimental density profiles for h, ~ 5 cm and h. = 10 cm: line with 
open circles, initial profile: line with closed circles, final profile (a) v = h. St 
= 205 mjp; (b) t = b + 120/i, St - 185 min. S indicates elevation of ob¬ 
served (dyed) interface 


Fig. 5. Schematic drawing of internal wave gage 

Fig 6. Theoretical and experimental waves at x = 6 = 61 cm 
with h = 15 cm and e = 0.05: solid line, theory: dashed line, 
experiment, (a) Free surface; (b) interface: /i, = 5 cm. C„ = 1.00 
cm. r, = 0.13. r, - 0.0(4. (c) interface: h, = 7.5 cm. Cn = - I 00 
cm. r, = 0.25. t, = 0.027; (d) interface: h, - 10 cm. Co = 1.08 
cm, r, = 0.12. r, = 0.012. 

Fig 7. Theoretical and experimental waves at i = h = 61 cm 
with h - 15 cm and e = 0.05: solid line, theory, dashed line, 
experiment, (a) Free surface; (b) interface It, - 5 cm. Cn = -LOO 
cm, t, = 5.13. r 2 = 0.53, (c) interface, h, = 7.5 cm, Co = 0.97 

cm. r, = 4.53, Tj = 0.50; (d) interface: It, = 10 cm, Co = 2.00 
cm. r, = 4 87, r, = 0.50. 
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Fig. K. Theoretical and experimental waves at x - b + 20 ; i: 

It = 15 cm. ft, = 5 cm. h = 61 cm, c — 0.05, U = 1.5 cm. 
r, = 0 12. t, = 0.012; solid line, free surface: dashed line, inter¬ 
face tat Linear theory: (b) experiment. 

Fit. 9 Variation of with r, and with t_. at i = b = 61 cm for It - 15 cm. 

t - 0 05. Solid curve represents surface wave calculation (for all /ip/n using linear theory Vertical 
slash indicates experimental data for surface waves 

FtO. 10. Theoretical variation of (.V,). na ,/£„ with r, and l(<V}) milv /(„)(/iiKi,| with t ; . at t = h. Dashed 
curve indicates asvmntotic results for h li t> 1 
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1. Introduction 


Studies of tsunami generation often employ simple models of the sea 
floor dislocations to permit tractable analytical solutions. Although 
these solutions provide basic insight into the generation process, they 
are incapable of producing explicit results for prototype events where 
both the spatial and temporal distributions of the sea floor dislocation 
may be quite complicated. Herein we exploit the apparent linearity of 
the generation process and demonstrate both the use and validity of the 
superposition principle to construct solutions for complex bed motions. 
Analytical and experimental results are presented for a monopolar dis¬ 
location (block upthrust or downthrow) with a complex time-displacement 
history. The time history used in the computations is obtained from an 
integrated accelerogram recorded at Pacoima Dam, near Los Angeles, 
during the earthquake of February 9, 1971. A complex spatial deforma¬ 
tion is not used in order to enable experimental verification of the 
analytical results. This is unfortunate since it appears that the 
details of the time-displacement history are not important for proto¬ 
type phenomena where the motion may be considered instantaneous. How¬ 
ever, it is important to note that the analysis treats both space and 
time variations in an identical manner; hence, confirmation of this 
approach for complex time variations strongly suggests analogous behav¬ 
ior for complex spatial variations. 

Finally, we examine and compare several alternative time-displacement 
histories for the mean motion. It is shown that the results for each mean 
motion can be unified by introducing a velocity as a descriptive param¬ 
eter which is based on the kinetic energy input of the moving bottom to 
the overlying fluid. 


2. Time-Displacement History 

To define a complicated time-displacement history with relevance to 
earthquake-induced ground motions, we have chosen an accelerogram for 
the vertical component of motion recorded at Pacoima Dam (near Los 
Angeles) during the earthquake of February 9, 1971. The accelerogram 
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component of ground acceleration recorded 2/9/71 0600 PST Pacoima Dam 
a, with integrated velocity and ground displacement. 













shown in Figure 1 may be integrated numerically using appropriate pre¬ 
cautions (e.g. , see Nigam and Jennings, 1968), to yield a time history 
of the vertical ground displacement. Results of the integration pro¬ 
cess for the "smoothed" velocity history and the corresponding ground 
displacement history are shown below the accelerogram in Figure 1. 

An analytical representation of the algebraically complex time- 
displacement history of Figure 1 over a finite record interval, 

0 < t £ T, may be constructed using a Fourier series of the form: 

N 

eft) = C + l C sin (K t + 3 ). (1) 

o , n n n 

n=l 

In (1) C n are the amplitudes of the Fourier components with C Q repre¬ 
senting the mean (permanent) ground displacment, K n are the component 
wave frequencies, and 3^ are the component phase angles. The accuracy 
of (1) in representing the integrated displacement of Figure 1 is 
determined by the number of components N retained in the Fourier sum. 
Results of computations with N = 18 are illustrated in Figure 2; this 
truncated sum will be adopted in the subsequent analysis. 

It should be emphasized that both the instrument characteristics 
and the numerical integration techniques used to obtain ground dis¬ 
placements from acceleration measurements necessarily distort (filter) 
information in long period components. In particular, the mean (and 
permanent) displacement of the integrated motion shown in Figure 1 (and, 
of course, its Fourier representation in Figure 2) is not expected to 
accurately model the actual permanent deformation. To compensate for 
this distortion, we may again exploit the superposition principle and 
add a nonsinusoidal component to the Fourier series representation of 
(1). As an example of this approach, consider the ramp motion in time 
of the form: 


C r (t) = C 0 t/T (2) 

during the time interval 0 < t £ T; other choices for the mean motion 
are discussed and compared in section 5. In practice, the sum of the 
mean component resulting from the Fourier synthesis, C 0 , and the addi¬ 
tional component Q 0 from (2) should be chosen to equal the actual per¬ 
manent ground offset. Hence, a general representation for a complicated 
time-displacement history becomes: 


C(t) 


c + 

o 


C 0 t/T 


• £ . 

n=l 


C sin 
n 


(K t 
n 




(3) 


3. Solutions of the Water Wave Problem 


Consider a two-dimensional (x,y) and incompressible ocean of uniform 
depth h initially in equilibrium with the earth's gravitational field g 
which acts in the negative y direction. At time t = 0 a section of the 
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sea floor begins to deform vertically with a time and spatial variation 
given by £(x,t). We seek the inviscid, irrotational, and barotropic 
deviations ri(x,t) of the ocean free surface from its equilibrium posi¬ 
tion. With the coordinate system at the equilibrium position of the 
free surface, the linearized description of motion in terms of a veloc¬ 
ity potential <p = 4>(x,y,t) is: 


4> xx (x,y,t) + 4> yy (x,y,t) = o 

(4) 

4> y (x,-h,t) = C t (x,t) 

(5) 

4> y (x,0,t) = n t (x,t) 

(6) 

4> t (x,0,t) = -gn(x,t) 

(?) 

where subscripted variables indicate partial differentiation, 
convenient to eliminate n(x,t) in (6) and (7) by combining to 
single condition for the velocity potential: 

It is 
yield a 

4> tt (x,0,t) + g<J) y (x,0,t) = 0 

(8) 

Using the Laplace transform in t and the Fourier transform in 
tions (4), (5) and (8) become: 

x, equa- 

? yy (k,y,s) - k^(k,y,s) = 0 

(9) 

^ y (k ,-h, s) = sc; (k , s) 

(10) 

^ y (k,0,s) + <jT(k,0,s) = 0 

(11) 

where the overbar of a function f(x,t) indicates: 


f(k,s) = | dx |°° e^ X e St f(x,t)dt. 

(12) 


-"O -oo 


Solving (9), (10), and (11) for <J>(k,y,sl and noting from (7) that 

n(k,s) = - (s/g)^(k,0,s) (13) 

we find: 

t 

r)(k,s) = s 2 C(k,s)/(s 2 + u) 2 ) cosh kh (14) 


where w 2 = gk tanh kh. Inverting the Laplace and Fourier transforms 
yields: 


n(x,t) 


r co r . 

_L f / lim 
2it J 

_oo 


1 

2fri 


|U+ir 

M-ir 


s 2 e~ lkx e st qk,s) ds | 
(s 2 + a) 2 ) cosh kh i 


dk 


(15) 
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Again, the final integration over wavenumber k must be evaluated numeri¬ 
cally. 


4, Comparison of Theory and Experiment 

As noted earlier, the simple block deformation of the sea floor was 
chosen to enable experimental verification of the analytical model 
developed for complicated time-displacement histories. The experimental 
facility used in these tests has been described in detail by Hammack 
(1972, 1973). Basically, the wave-maker consists of a rectangular 
piston in the bottom of a wave tank (and spanning its width) whose 
motion is controlled by an electro-hydraulic-servo system. The servo 
system converts a time-voltage command signal into a proportional verti¬ 
cal displacement of the piston. For the experiments reported herein, 
the piston length in the direction of wave motion is b = 61 cm while the 
quiescent water depth above the piston is h = 10 cm. Before presenting 
results of the tests, we describe the motivation for choosing other 
experimental scales. 

A "global" time scale for the forcing of the overlying ocean by the 
sea floor is the period T. The appropriate time scale for the gravita¬ 
tional response of the long barotropic wave modes is b/fgh)^ which cor¬ 
responds to the time required for waves to escape the generation 
region. For prototypical earthquakes the ratio of the forcing and 
response time scales, t = Tfgh^/b, termed the time-size ratio, is small 
so that details of the time-displacement history generally are not 
important. However, our interests herein require that the details of 
the temporal motion have a significant impact on the generated wave 
structure. Hence, the period T for the experimental tests must be 
scaled so that t exceeds unity; in fact, for the experiments a period 
T = 4 secs was chosen which yields t = 6.5. Previous experiments by 
Hammack (1973) also indicate that the generation process for proto¬ 
typical tsunamis is linear and that nonlinearity remains insignificant 
for vertical displacements which do not exceed about 20% of the over- 
lying ocean depth. This criterion is adhered to in the experiments by 
restricting the instantaneous displacement of the piston to less than 
2 cm. 

In the first test we examine experimental and theoretical results 
for the time-displacement history shown in Figure 2 using experimental 
parameters T = 4 secs and a chosen permanent (mean) displacement of 
C Q = 0.38 cm, with the amplitudes and frequencies Kj, of the eighteen 
Fourier components used in Figure 2 scaled appropriately. (In Figure 2 
the corresponding parameters of the actual ground displacement are 
T = 40 secs and C Q = 7.63 cm.) The scaled Fourier components are then 
summed and the result is converted to an analogue (time-voltage) 
signal which is used to command the wavemaker. The resulting wave 
motion at the leading edge of the piston (x = b) is measured, and the 
results are shown in Figure 3. Theoretical results at x = b are 
evaluated for each of the eighteen (scaled) Fourier components accord¬ 
ing to (17) and summed to yield 









the results also are shown in Figure 3. The excellent agreement between 
the predicted and measured data is self-evident. 

Similar results are shown in Figure 4 where a ramp mean motion with 
an amplitude of Co = 1.33 cm is added to the Fourier synthesis of Figure 
2. The total permanent displacement of the piston is Co’ = + Co = 

1.71 cm which has been used to normalize the measured and theoretical 
wave amplitudes. The theoretical result is equivalent to (25) with an 
added component for the ramp computed from (24). Again, the agreement 
between measured and computed data is excellent with the wave structure 
clearly showing the added volume (mass) resulting from the enhanced mean 
displacement. 


5. A Comment on Mean Motions 

In previous studies (Hammack, 1972, 1973) two additional models for 
the mean displacement of a block section of the sea floor have been 
examined. These time histories are: 

a. exponential: C e (t) = C Q [1 - exp (-1.1t/T e )] 

b. half-sine: C g (t) = C Q [(1 - cos irt/T s )H(T s -t)/2 + H(t-T s )] 

and we repeat for completeness the mean motion introduced here: 

c. ramp: C/t) = C Q [tH(T r -t)/T r - H(t-T r )]. 

The three mean motions listed above span a wide range of displacement 
characteristics. We note that the choice of characteristic time scales 
T e , T s and T are, in fact, arbitrary to a certain extent even though 
"natural" choices are apparent. (This flexibility is most obvious for 
the exponential motion where T e was chosen for experimental convenience 
to represent the time for two-thirds of the displacement to occur.) It 
has been found that the properties of waves generated by these motions 
correlated strongly with the time-size ratio x based on these time 
scales. For example, the maximum wave amplitude, say q 0 , occurring at 
x = b when normalized by the permanent displacement £ 0 exhibits a simple 
and similar functional dependence on x for each bed motion. With the 
size scale b/h of the dislocation fixed and for x<<1, the normalized 
amplitude, r\ 0 /C, 0 , teaches a maximum value of one-half for all size scales 
exceeding unity. Bed motions with x<<1 are termed impulsive. For x 
very large, termed creeping generation, Ho/^o decreases at a rate which 
is inversely proportional to x. The constant of proportionality for 
creeping generation varies with the specific choice of the characteris¬ 
tic time scale for the mean motion. Since it is unlikely that any of 
these mean motion models is "correct" from a geophysical point of view, 
there is a need to seek a unification of results by generalizing the 
concept of characteristic time scale. One generalization which closely 





where the complex inversion integral for the Laplace transform has been 
used. In (15) ui is the wave frequency (w = is) and k is the wavenumber. 
Explicit results for specific deformations of the sea floor £(x,t) will 
be developed now. 

3.1 Solution for a single Fourier component 

Consider a block section of the sea floor of length 2b whose time- 
displacement history corresponds to a single Fourier component of (1). 
With the coordinate system centered above the block section, we have 

C n (x,t) = C n H(b 2 - x 2 )[sin(KT+6 n )H(T-t) + sinf^t + 0 n )H(t - T) ] (16) 

where H( ) is the Heaviside step function. Finding the transform of 
(16), substituting into (15) , performing the integration around the 
Bromwich contour, taking only the real part of the resulting integral, 
and noting that the integrand is an even function of k, we find: 

%<*■« ■ ^ f ° ° f l l rh 1 fa r b {**B-H(t-T)[C.D.E]} dk (17) 

o 

where 


A = 

sin B n [(to 2 sin ait - 

K 2 sin 
n 

K n t)/(W 2 - K n 2 )] 

(18) 

B = 

K cosg [(a) sin out - 
n n 

k n sin 

K n t)/(0J 2 - K n 2 )] 

(19) 

C = 

sln(KT . 6 n )j 

[[u 2 cos to(t-T) - K n 2 cos K n (t-T)]/(U) 2 -K n 2 ) 

} (20) 

D = 

K n cos(lt „™n>i 

r 

[gj sin 

oj(t-T) 

- K sin K (t-T) ] / (uj 2 - K 2 ) 1 
n n J J v n J 

(21) 

E = 

- sin (K T + 8 ) cos 
n n' 

u(t-T) 

. 

(22) 


The final integration over wavenumber k in (17) is obtained by numerical 
quadratures. 

3.2 Solution for ramp mean motion 

The ramp time-displacement history of (2) for the block deformation 
is described by 

C(x,t) = H(b 2 - x 2 ) [tH(T - t)/T ♦ H(t - T)]. (23) 

Following the same procedure outline in section 3.1, we obtain (again 
after considerable algebra) the water surface motion n r due to the ramp 

V-o ■ ~ir Prlos-ra- (Jrji 51 " “< t - T » dk -' 24 ) 
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RESULTANT WATER WAVE AT x = 



Figure 3. Theoretical (-) and experimental (—) wave records at 












Theoretical (-) and experimental (-—) wave records at x = b with ramp mean motion. 










Theoretical and experimental variati 













produces the desired unification is the following. Consider a velocity 
V, where V 2 is the kinetic energy imparted to the fluid by the sea floor 
divided by one-half the total mass of fluid displaced during bed motion 
PspbC 0 ); then V is given by: 

1 

v = ~ J C t 2 dC (26) 

and a corresponding time scale T = C 0 /V. In terms of this time scale 
the time-size ratio becomes 

T* = T(gh) % /b = C 0 (gh) i5 /bV (27) 

The variation of p 0 /C 0 at x = b for b/h = 6.1 with t* is shown in Figure 
S for each of the mean motions; both theoretical and experimental 
results are presented. For all of the experimental data we have taken 
£ 0 /h < 0.2 to avoid significant nonlinear effects. The collapse of 
results for such a wide range of mean motion characteristics shown in 
Figure 5 is good although a small spread still exists. 


6. Conclusions 


We have demonstrated both the application and validity of a strategy 
which employs multiple uses of the superposition principle to develop 
theoretical solutions for waves generated by sea floor motions with com¬ 
plicated time-displacement histories. Although a more useful test for 
prototypical phenomena would utilize complicated spatial distributions 
for the sea floor dislocation, the tests herein were restricted to 
simple block dislocations due to experimental limitations. However, it 
is emphasized that the solution method does not distinguish between 
space and time, and the validity established herein strongly suggests 
that the methods could be extended to complicated spatial deformations. 
Finally, we have demonstrated that wave properties (in particular the 
maximum amplitude of waves escaping the generation region) for a wide 
range of mean motion characteristics may be (almost) collapsed into a 
single functional relationship in terms of a time-size ratio based on an 
average vertical velocity of the sea floor obtained from energy consid¬ 
erations. 
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SMALL-SCALE OCEAN WAVES 

J. L. Hammack 

Department of Civil Engineering, University of California - Berkeley, California 

I. — Introduction. 

Presented below are three lectures on various aspects and types of ocean 
waves. In all cases we are concerned with small-scale waves for which gravi¬ 
tation provides the dominant restoring force. The relevant length and time 
scales of these waves are small compared to the Kelvin-Fossby radius of defor¬ 
mation and inertial period, respectively, and Coriolis effects may be 
neglected. In addition, the small length scales (relative to Earth's radius) 
permit us to ignore the curvature of the ocean surface and adopt a plane- 
earth approximation. Various other approximations will be adopted as 
appropriate in order to obtain model equations which are tractable analytically. 
(Fortunately, these tractable models appear to remain relevant for geophysical 
phenomena!) In particular, we are interested in weakly nonlinear systems in 
which the nonlinear effects manifest slowly; both nonlinear self-interacting 
and resonant-triad systems are considered. 

An outline of these lectures is as follows. In Part II we examine two 
nonlinear model equations for the evolution of gravity waves: the Korteweg- 
deVries (KdV) equation for long barotropic and baroclinic waves and the 
nonlinear Schrodinger (NLS) equation for short barotropic waves. After a 
brief review of the asymptotic (large time) solution of these equations by 
inverse scattering theory, experimental data are presented which demonstrate 
the reality of soliton predictions—at least on laboratory scales. In Part III 



2 


we exploit the exact solution of the KdV equation and its linear approximations 
in order to derive rather precise criteria for modelling the evolution of 
long-wave initial data on geophysical or laboratory scales. The modelling 
criteria of Part II are then applied to a typical oceanic tsunami in order 
to choose relevant model equations for propagating the (barotropic) wave from 
its generation region to specific target sites. Finally, in Part IV we 
examine the nonlinear excitation of "edge" waves nearshore by linear wave 
trains incident from deep water. Both theoretical and experimental results 
are presented which document the excitation of two progressive edge-wave inodes 
through a nonlinear resonant-triad interaction with the wave reflected after 
normal incidence from offshore. 

Much of the material in Parts II and III has appeared previously in 
the literature. Hence, many of the details concerning experimental equipment, 
procedures, and analysis is omitted; a thorough discussion may be found in 
the cited references. The material presented in Part IV on edge waves is 
recent and not presently available in the literature. For clarity, a more 
detailed discussion of this material is presented. We also note here that 
the notation between Parts II, III, and IV is not necessarily consistent. 

II. —Water Wave Solitons. 

1. — Long Gravity Waves 

Consider two fluid layers with uniform densities resting in a 

gravitationally stable configuration (see Figure 1) on a horizontal and 
impermeable bed of infinite lateral extent. The upper layer possesses a 
free surface along which pressure is constant; surface energy effects on 
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and the .fluid-fluid interface are negligible. Required are the two- 

dimensional, inviscid, irrotational motions which are bounded everywhere and 

evolve from given initial distribution of velocities, free surface deformation, 

and interfacial deformation—the classical water wave problem. In particular, 

we are interested in the barotropic displacement of from its static 

equilibrium position and the baroclinic displacement r\ of from its 

equilibrium position. To derive the KdV equation as the appropriate model 

for both displacements from the governing equations, the following assumptions 

are necessary. First, characteristic wave lengths k ^ (an inverse wave 

2 

number) must be long relative to the total fluid depth; i.e., (kh) << 1, 

so that dispersive effects are weak. Second, characteristic wave amplitudes 

a are small relative to the total depth; i.e., a/h << 1, so that nonlinearity 

is weak. Third, both weak effects of dispersion and nonlinearity are approxi- 

2 

mately equal; i.e., e = a/h 'v (kh) . For simplicity, we will further assume 

at the outset that density differences in our two-layer ocean are small (as in 

its geophysical counterpart) so that A = (p^ - << 1- Expanding the 

dependent parameters in the governing equations in a power series in e and 

introducing the multiple time scales tg * t, t =» et, . . ., one finds at 

leading order that the evolution equations for the initial wave are hyperbolic 

(nondispersive) and linear. At this order, which corresponds to the fast time 

scale tg, an initial disturbance decomposes into four modes consisting of 

left- and right-running, barotropic and baroclinic modes. The phase speeds 

C, and C of these modes are 
1 2 







2 

(la) Barotropic (surface) waves: = gh 

2 

(lb) Baroclinic (interfacial) waves: = g'h^h^/h 


where we have invoked the Boussinesq limit A -*■ 0 with gA = g* << g remaining 

finite. All wave modes propagate with permanent form and do not interact 

with each other or themselves; the baroclinic inodes propagate much slower, 

1/2 

0(e ), than the barotropic modes. 

2 

At the next order (s ) weak nonlinear effects and dispersion occur. 

Each wave mode experiences a self-interaction on the slow time scale t = et, 
but no interactions between modes occur due to their rapid separation by the 
phase speed differences of (1). The self-interaction of the right-running 
barotropic mode is governed by a dimensional equation of the form 


( 2 ) 


n i + c i n i + + h h \ 

t X X XXX 


= 0 


A more convenient choice of nondimensional variables for describing these 


waves are 


<3> < 


X = (x - C t)/h, r = 7 (g/h) 1 ' /2 t 

J- 6 


f(x»T) = jn(x,t)/h 


V 

In terms of these variables, (2) reduces to the KdV equation with the common 
form 


(4) f + 6ff + f = 0 

t X XXX 


In a similar manner, the dimensional equation governing the self- 

i 


interaction of the right-running baroclinic waves is 







n_ + c n + -c (--— )n_n~ + tCJi h n 

2^ 22 22 h_ h, 22 62122 

t x 2 1 x xxx 


Introducing the normalized variables 


x = (x - c 2 t)/(h 1 h 2 )‘ , ‘ 

1 1/2 

(6) < t = -Hg'/hJ^t 
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, 3 1 1 . 

f = -(--7—)n 

2 h 2 hj 

equation (5) reduces to the KdV equation given by (4). 

The asymptotic solution of the KdV equation for arbitrary initial data 
f(X/0) H f Q by inverse scattering theory has been described in detail by 
Segur'*’ inter alios . Here, we briefly list features of the asymptotic solution 
to be illustrated in the laboratory experiments. 

i) An initial disturbance evolves into a finite number of permanent 

waves (solitons) ordered by their amplitude. When the solitons are 
well separated, the local shape of each is given by 

2 1/2 

(7) f = asech {(a/2) (X “ X Q - 2ax)} 

where a and Xq are constants. The rank ordered solitons are 
followed by a dispersive train (radiation) of oscillatory waves, 
ii) The number N of solitons evolving from initial data of finite 

extent, say f = 0 for x < X^ and x > X 2 ' i- s equivalent to the number 
of zeros of $ for X > X^ where <J> satisfies 


±A + 
dx 2 


2 + V XH = ° 


l( xj» lM (xJ = 0. 

L dv 1 





iii) When the net volume V (or mass) in the initial wave is finite and 


positive, i.e. 


(9) 


V = 


f Q (x)dx > 0 


at least one soliton emerges. 

iv) When f < 0 everywhere, no solitons emerge and the asymptotic 
solution consists (only) of the radiation components, 

v) Two other important classes of data are those for which V < 0 with 
f q 4^ 0 for all x and those for which V = 0. No general statements 
regarding the asymptotic solution for these cases is provided 
theoretically. Experiments demonstrate that the evolution of 
solitons depends on the detailed structure of the initial data. 


1.1 Experimental Procedures In order to illustrate the applicability 

cf the KdV equation as a model for long water waves, a series of experiments 

is conducted in a laboratory wave tank 31.6 m long, 61 cm deep, and 39.4 cm 

wide. For studying the barotropic motions in a system with small density 

differences, the density stratification plays no role; hence, in these 

experiments a uniform density fluid of depth h is used. A detailed descrip- 

2 

tion of these experiments is given by Hammack and Segur. In order to study 
the evolution of baroclinic waves, the tank is stratified with fresh water 
and brine. Details of these experiments can be found in Hammack^ as well as 

4 

Hammack and Segur. 

In the long wave experiments reported here, the wavemaker consists of 
a rectangular piston located in the tank bottom adjacent to an upstream 
end-wall. The piston spans the tank width and has a length b in the direction 






of wave motion (see Figure 2); lengths of b = 30.5 cm and 61 cm are used. 

The vertical motion of the piston is controlled by am electro-hydraulic-servo 
system and completely user specified. In a typical experiment, the piston 
is moved for a finite time interval. Differences between the initial and 
final position of the piston permit the net volume V in the generated wave 
train to be calculated. Various initial data are generated by varying the 
time-displacement history of the piston. 

Both surface and internal wave amplitudes are measured at fixed loca¬ 
tions (stations) along the tank. We note here that differences between the 
temporal variation in wave amplitude at a fixed spatial position and the 
spatial variation and wave amplitude at a corresponding fixed time are small, 
0(e), and neglected in all calculations. 

1.2 Results for Barotropic Wave Evolution 

1.2.1 Initial data with V > 0. 

Figures 3 and 4 illustrate the evolution of two barotropic waves with 
a net positive volume. Normalized wave amplitudes are presented at four 
succeeding stations along the tank in a coordinate system which moves with 
the linear (nondispersive) speed C^. Note that the leading portion of the 
wave system appears at the left in these figures. Shifts of the waves to the 
left (right) at succeeding stations indicate phase speeds greater (less) than 
in this coordinate system. 

The initial wave at (x - b)/h = 0 in Figure 3 is rectangular and positive 
and appears to sort itself into three positive waves identified by separate 
crests (local maxima) during propagation. These three waves are rank ordered 
by amplitude and are followed by a weak train of dispersive waves whose speed 
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is much less than C^. In order to examine the local shape of the leading two 
waves, theoretical soliton profiles defined by (7) have been superposed on 
the measured data at the last measurement station (using the measured wave 
amplitudes for specifying the parameter a). These profiles are not extended 
into regions where the two waves are still interacting strongly with adjacent 
wave structure. Clearly, the measured waves appear to be locally KdV solitons. 
Further evidence that the three lead waves at the last measurement station 
are solitons is provided by numerically integrating (8) using the wave profile 
at (x - b)/h = 0 as the initial data f^. Computations indicate that N = 3 
solitons should evolve in agreement with the observed pattern of evolution. 

The evolution of a more complicated initial wave with V > 0 is shown 
in Figure 4 where a leading negative wave is followed by a (larger) positive 
wave. After only twenty depths of propagation, the positive wave has 
separated into three separate crests while the negative wave remains essen¬ 
tially unchanged. During subsequent propagation the three labelled crests 
of the positive wave appear to retain their integrity as they progress thro’ 
the leading negative wave and emerge at the front of the wave train. At the 
last measurement station, labelled waves 1 and 2 clearly resemble KdV 
solitons. The third wave is still interacting with the once-leading negative 
wave at the last station and cannot be unmistakably identified as a soliton. 
However, computations with (8) using the wave at (x - b)/h = 0 as f^ yield 
N = 3 strongly suggesting the third wave is indeed a soliton. We note that 
other experiments also indicate that solitons evolve from the positive waves 
in the initial wave and can be identified long before asymptotic conditions 


fmmm - 







are achieved. 


1.2.2 Initial waves with V < 0 


Results for the evolution of two barotropic wave systems with V < 0 are 
shown in Figures 5 and 6. The initial wave in Figure 5 is the negative counter¬ 
part of the experiment shown in Figure 3 where three solitons appeared to 
evolve. Over similar distances of propagation, no solitons appear to evolve 
in Figure 5 — just as expected for initial data with f < 0 for all x- 
Instead, a negative wave evolves whose frontal slope decreases and lengthens 
with time. This lead wave is followed by a train of strongly dispersive 
waves with phase speeds much less than C^. In fact, the wave structure of 
Figure 5 represents the radiation solution of the KdV equation as shown with 
quantitative tests by Hammack and Segur. 5 

Further evidence for the evolution of complicated waves with V < 0 is 
shown in Figure 6. In this case the lead positive wave evolves into one or 
possibly two waves while the negative wave evolves in a manner similar to 
that observed in Figure 5. Computations using the initial wave in (8) indicate 
that one soliton should develop asymptotically. Hence, the second labeled 
wave is expected to eventually disappear as it appears to be doing. 

1.2.3 Initial waves with V = 0 

Figure 7 illustrates the evolution of initial data with V = 0 consisting 
of a positive wave preceded and followed by negative waves. Applications of 
(8) to the initial wave suggest that two solitons should evolve. The large 
positive wave in the initial data quickly separates into two crests which 
appear to migrate through the lead negative wave during subsequent propagation. 
At the last station, labeled wave 1 has progressed to the front and clearly 
has the shape of a soliton everywhere. The second labeled wave still appears 
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to be interacting with the frontal slope of the once leading negative wave. 

A trailing grain of dispersive waves similar to that of Figure 5 also evolves. 

1.3 Results for Baroclinic Wave Evolution Ihe evolution of baroclinic 

long waves at the interface of a stratified fluid with h^ = 45 cm, h^ = 5 cm, 

and A = 0.05 is shown in Figure 8. Wave amplitudes recorded at seven 

stations along the tank in a coordinate system that moves with the linear 

speed are presented. The initial wave is positive, and according to 

calculations by (8), should evolve two solitons. Ihe observed evolution is 

in agreement with this prediction. (A strongly damped train of oscillatory 

waves has been omitted in Figure 8.) A more quantitative comparison of the 

lead wave profile with the theoretical shape (7) at the last four stations of 

measurements is shown in Figure 9. The agreement with the measured data in 

Figure 9 is excellent. Further results for baroclinic wave evolution are 

4 

presented in Hammack and Segur. 

2. — Short Gravity Waves 

We now turn our attention to the evolution of short gravity waves; i.e., 
waves whose characteristic length k ^ is comparable to or greater than the 
local ocean depth; only barotropic modes are considered. In fact, our 
interest in solitons requires kh > 1.36 which we adopt as the definition of 
the short-wave regime herein. In order to derive the nonlinear Schrodinger 
(NLS) equation as an evolution model, we consider a wave train with a 
dominant and identifiable (mean) wave number k^. As for long waves, we require 
weak nonlinearity, which is now characterized by the wave steepness e = ak Q << 1. 
The wave system is permitted to have weak modulations such that the variation 
in wave number 6k is small; i.e., 5kA 0 <<: 1* In other words,, we are concerned 






with "narrow-band" wave systems that are weakly nonlinear. Thirdly, we 

postulate a balance of both small effects so that e = ak^ *\» 5k/k g . Multiple 

scale analysis of the governing equations again yields a heirarchy of 

problems at different orders of e. At lowest order we recover the linear 

£ 

dispersive waves of Stokes with an amplitude a that is constant. At the 

next time scale, t = st, the wave amplitude is modulated, and we find that 

the amplitude modulations propagate with the linear group speed C . 

2 

Continuing to the third order with a time scale t^ = e t, one finds that the 
complex amplitude modulation a must satisfy the NLS equation. Defining non- 
dimensional coordinates as 

f 

X = ek 0 ( x - C g t) . 

(10) / t = e 2 (gk 0 ) 1/2 t 

A = k 0 2 (gk 0 r 1/2 a 

V 

the NLS equation takes the form 

(11) iA + AA +v|a| A = 0 

XX 

where A and v are known functions of the water depth h, gravitation g, and 
carrier wave number k Q (or frequency u) Q ). Details of the derivation for 
finite depth can be found in Hasimoto and Ono^ inter alios . 

Like the KdV equation, the NLS equation cam be solved exactly for 
arbitrary initial data by inverse scattering theory. When k Q h > 1.36, 


envelope solitons can occur; the one-soliton solution of (11) is 
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2 \ ' 2 
(12) A = a — sech{ax}exp{iAa t} 

where a is an arbitrary constant related to the envelope amplitude. In general, 
initial data of finite extent will evolve a finite number of envelope solitons, 
rank ordered by the group velocities of their dominant carrier waves and 
embedded in a dispersive train (radiation) of oscillatory waves which decays 
in amplitude with time. It is important to note that the speeds of both the 
solitons and the radiation components are not impacted by nonlinearity at this 
order. Hence, unlike their long wave counterparts of the KdV equation, these 
solitons and radiation components do not separate with time. This suggests 
that in order to observe clearly the evolution of envelope solitons even in 
a contrived laboratory experiment, one must design the initial data such that 
the wave content of the radiation spectrum at the dominant frequency of the 
soliton carrier waves is small. Alternatively, one must observe evolution 
until the inviscid decay of the radiation (by frequency dispersion) combined 
with viscous decay has progressed sufficiently. 

Although a quantitative comparison of theory and experiment for the 

NLS equation analygous to that of the KdV equation has not been performed, 

0 

qualitative tests by Yuen et al provide evidence for its applicability to 
narrow band systems. We present in Figure 10 an (unpublished) experiment 
which illustrates the evolution of an initial wave packet whose envelope 
amplitude and carrier wave frequency are modulated; the mean wave frequency 
is = 0.6 sec. The experiment is conducted in a tank approximately 
50 m long, 2.5 m wide, and a water depth of h = 1 m. After 45 m of propaga¬ 
tion, the initial packet evolves into a collection of wave groups which are 
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ordered by the group speed of their dominant carrier wave. Energy is concen¬ 
trated in the group with carrier waves at the dominant period u> 1 * 0.6 sec of 
the initial data. The envelope of this group clearly has the shape of the 
soliton profile (12) shown superposed on the measured wave. The envelopes 
for the other wave groups do not agree with (12) , and apparently these groups 
represent the radiation components in the initial data. 

Although testing of the NLS equation is not complete, there is growing 

9 

evidence (see Yuen and Lake ) that the dynamics of short-water waves is much 
more complicated than that of long waves,and more complicated models are 
required. The restriction of the NLS equation to narrow band systems and its 
probable) instability to transverse perturbations appear especially severe for 
many geophysical applications. 

XXI. — Modelling Criteria for Long Water Waves. 

Even with all the assumptions implicit in the classical water-wave 
problem introduced in Part II, the general equations remain intractable analyt¬ 
ically. Hence, we resort to further approximations such as those required to 
yield the KdV and NLS equations. Generally, approximations are formalized by 
perturbation expansions in terms of a small parameter(s). In application of 
these approximate models, questions naturally arise as to when the inequalities 
used in ordering the physics are actually satisfied; i.e., how small is small? 
For evolution models, the relevant question distills into "during what time 
interval does a particular approximation correctly represent the general 
solution?" Closely related to this topic are questions of similitude between 


laboratory models and their geophysical-scale counterparts. The answers to all 
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of these questions involve determination of nondimensional parameters which 
characterize the phenomenon in question. In practice these parameters are 
often deduced through the ingenious use of dimensional analysis, by examina¬ 
tion of the governing equations, or using the most preferred but least available 
method, examination of exact solutions of the general and approximate model 
equations. Here we examine the KdV equation and exploit exact solutions of its 
various approximations in order to develop rather precise criteria for 
modelling the propagation of long water waves. These criteria are then applied 
to the problem of tsunami propagation across ocean basins—the topic which 


precipitated the analysis. Further details are presented by Hammack and Segur.' 
1. — Analysis. 


In order to obtain definitive results, we adopt the following point of 

view. If the initial wave f is sufficiently smooth and localized, and if 

2 

= a/h and = (kh) based on the initial data are both small, then the 
KdV equation (4) will be the approximate model eventually. (In fact, we know 
that (4) is appropriate when e = 0(e ^).) If and eare not the same 

order of magnitude initially, then simpler forms of the KdV equation may be 
applicable for some time interval, e.g.. 


f = 0 (linear, nondispersive model) 


f + f 
I XXX 

f + ff 
I X 


0 (linear, dispersive model) 

0 (nonlinear, nondispersive model) 


Ursell showed that the nondimensional parameter 


U = = ak 2 /h 3 
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which we shall call the Ursell number, is an important indicator to determine 
which of these models is appropriate. One should use (14) , (4), or (15), 
depending on whether U << 1, U = 0(1), or U >> 1, respectively. Ursell also 
provided evidence that U is time dependent and will tend to an order-one 
limit so that all waves in this category eventually propagate according to 
(4) as postulated previously. Although there is general agreement with 
Ursell's results, there is disagreement on how to interpret "order unity" and 
how to define the relevant length scale k ^ for the evolving waves. To make 
his criteria more precise, we examine here initial data for which << << 1 

initially so that the linearized models (13) and (14) are the relevant approx¬ 
imations of (4). The dimensional length scale k Q 1 for the initial data is 
defined as its overall length, and the dimensional amplitude scale is a Q . 
Further, we will focus on criteria for modelling the leading wave only. 

Results for the trailing wave structure and the nonlinear approximation (15) 
are given in Hammack and Segur. 10 

1.1 Linear Dispersive Theory Since nonlinearity is small for the 
postulated initial data, we assume a small parameter y « 1 exists and seek a 
formal series solution of the KdV equation (4) in the form 

(17) f = yf L + u 2 f 2 + 0(y 3 ) 

(The required definition of y will come out of the analysis.) Substituting 
(17) into (4) yields a heirarchy of problems. 


(18) 

(f, ) + 

(f,) 

1 X 

1 XXX 

(19) 

(f,) + 

(f,) 

2 T 

2 XXX 


etc. 



= 0,yf 1 (x,0) = f Q 

= -6W ,f 2 (X.O) = 


0 
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The solution of (18), which is equivalent to (14), is well known: 

CO 

(20) Uf I = 27 { f 0 ( K ) ex P [*■ (<X + < 3 T)]dK 

—CO 

whe re 

OO 

(21) f Q (<) = { f Q (x)exp(-iKx)dx 

— OO 

is the Fourier transform of the initial data. For practical reasons, we are 
most interested in the asymptotic form (t -*• “>) of the solutions for both the 
linear and nonlinear models. As t •> » with (x|/t ->-0 in order to remain at the 
wave front, the asymptotic form of (20) is 


(22) Mf L (x»T) = f Q (0)(3T)‘ 1/3 Ai(C) - if^(0)(3T)" 2/3 Ai'(C) 

- |f^(0) (3T) _1 Ai"(S) + 0 [ (3 t) ~ 4/3 ] 

1/3 

where C = X/(3t) and Ai(£)is the Airy function. The coefficients in (22) 
have simple interpretations: 


(23a) 


(23b) 


(23c) 



where V is the nondimensional volume of the initial wave, is an Ursell 
number for the initial wave, and 8^, 8^ are constants that depend on the 
details of the initial wave. Thus (22) becomes 







(24) Uf 1 (X.T) = V(3T)~ 1/3 {Ai(C) + 6 1 (^r) (3T)" 1/3 Ai' (?) 

+ (3T)“ 2/3 Ai"(C) + . . .} 

From (24) it follows that the time required for this representation to become 
asymptotic (second term smaller than first, etc.) is at least 


(3t) 1/3 » U Q /|v| 


A particular solution of (19) is 


! f 2 (x»f) = - { / pf^U/Tjdz}^ 


and defining 


1 , x 

Q(x) « ~ {/ f 0 u)d Z r 


V — oo 


the homogeneous solution of (19) is 

X 2 » „ 

(28) U 2 f (X.t) = V 2 / Ai (z) dz + J~f (0(k) ~ (i*c> _1 > 


exp [i (icy + < t) ] dK 


Thus, as t + », jx|/ T 0 


U 2 f (X,t) = V 2 {- [ / Ai(z)dz] 2 + / Ai(z)dz + 


Comparing (24) and (29) , the appropriate definition of y is found to be simply 


V = V 











In other words, linear dispersive theory along with its large-time asymptotics 
requires the dimensionless wave volume V to be small. (Note that V is independent 
of time.) Since the solution of (29) remains 0(1) as t -*■ <*>, the series (17) 
cannot remain asymptotic after 

(31) (3t) 1/3 % |v| -1 

Thus, if the initial wave has a dimensionless volume V (assumed ^ 0) and 
Ursell number Uq based on its initial (given) dimensions, asymptotic linear 
dispersive theory is valid in an interval no longer than 

(32) Uq/M « (3t) 1/3 « l/|v|. 

1.2 Nondispersive Linear Theory. Since our initial wave is postulated 
to be linear, then U Q << 1 and necessarily [Vf «< 1, as well. The derivation 
of the KdV equation outlined in Part II indicated that linear nondispersive 
theory (13) occurs on the first time scale of evolution. Hence we may write 

(33) 3t « e 2 = (k Q h) 2 = (U Q /V) 2 

for (13) to be applicable. 

2. Summary and Application of Criteria to Tsunamic Propagation 

In summary, we have examined the evolution of long-wave initial data which 
is parameterized by an initial Ursell number U Q << 1 and volume (or mass) 

^ <<<: ma y model the evolution of this initial data using linear nondispersive 

theory (f =0) during a time interval 
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0 < 3t << Ti 


The next relevant model is linear dispersive theory (f + f =01 

t XXX ’ 


which is valid during 


Ti << 3x << T 2 = U 0 3 


with its asymptotics becoming valid during 


t 2 << 3 t << t^ = y2. 


Subsequently, we must use the KdV equation (f + 6ff + f =0) for 

t X XXX 


3x » T - 


(The question of the time scale for applicability of KdV asymptotics is discussed 

by Hammack and Segur^.) As an example of the application of these results, 

. . . 12 
consider the dimensional scales adopted by Carrier for major tsunamis which 

4-1 5 

impact entire ocean basins: a Q = 10 ft., h = 1.5 ' 10 ft., k = 2 * 10 ft. 

Then U Q t 0.1 and V % 0.01 and linear nondispersive theory is valid for times 

2 

corresponding to propagation distances (using = gh) of 


0 < x << 600 miles. 


Linear dispersive theory is appropriate for 


600 miles << x << 6,000 miles 


indicating that dispersion may affect this wave over much of typical ocean 
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trajectories. Even so, linear asymptotics do not apply until x >> 6,000 miles, 
which exceeds the length of realistic trajectories. Hence, linear dispersive 
asymptotics and the KdV equation are not required for describing the lead wave 
of the tsunami discussed here. (The KdV equation may be used, but it is 
unnecessarily complicated.) Of course, other factors may be required to 
accurately model long-term tsunami propagation (such as the variable bathymetry 
along trajectories and three-dimensional spreading of wave energy); the intent 
of the analysis here is only to develop insight into the relative importance 
of dispersion and nonlinearity. 

IV. — Excitation of Standing Edge Waves on Beaches 

1. — Introduction 

In recent years considerable attention has been focused on the occurrence 
of ocean waves which become trapped and capable of concentrating energy in 
localized regions (wave guides). Wave trapping can occur whenever gradients 
exist in a parameter which affects the wave's phase speed, e.g., Coriolis parameter, 
Brunt-Vaisala frequency, current speeds, and water depth. Stokes provided 
the first theoretical evidence of trapping for surface (gravity) waves near the 
shoreline of a plane sloping beach. Stokes found a normal mode solution for the 
(barotropic) departure n of the water surface from its static equilibrium position 
(see figure 11a) of the form 

(38) n (x,y,t) = a sin 6 exp (-k y cos 8) cos k x sin w t 

m m mm 

where -tan 8 is the beach slope, x points in the longshore direction, and 
y points offshore. The (linear) dispersion relation for the various longshore 
mode number m = 1,2, ... is 




- 19 


(34) 


0 < 3t << t i 


fn 1 


The next relevant model is linear dispersive theory (f + f = 0) 

t XXX ' 

which is valid during 


(35) 


T 1 << 3t << 7 2 = 


V 


with its asymptotics becoming valid during 


(36) 


T2 << 3x << T. 


V' 


Subsequently, we must use the KdV equation (f + 6ff + f =0) for 
" t x XXX 


(37) 


3x >> T. 


(The question of the time scale for applicability of KdV asymptotics is discussed 
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by Hammack and Segur .) As an example of the application of these results, 


12 

consider the dimensional scales adopted by Carrier for major tsunamis which 


4-1 5 

impact entire ocean basins: a Q = 10 ft., h = 1.5 ' 10 ft., k Q =2-10 ft. 


Then Uq 'v 0.1 and V ^ 0.01 and linear nondispersive theory is valid for times 


corresponding to propagation distances (using = gh) of 


0 < x << 600 miles. 


Linear dispersive theory is appropriate for 
600 miles << x << 6,000 miles 

indicating that dispersion may affect this wave over much of typical ocean 
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(39) to = gk sin B 

m m 

where k = imr/b is fixed by the beach width b. Note that the Stokes' mode 
m 

of (38) is periodic in the longshore direction with crests pointing offshore 

which decay in amplitude with an e-folding distance y = (k cos B) -- * - . 

e m 

The crest amplitude is maximum at the shoreline with a magnitude a sin B 

where a is referred to as the run-up amplitude. (Of course, run-up phenomena 

cannot be represented by linear solutions; however, the reality of run-up is 

acknowledged since a is the up-beach length of the horizontal projection of the 

14 

vertical shoreline amplitude.) Ursell demonstrated that the Stokes’ mode 
is only the lowest (j = 0) of a discrete set j = -0,1 ... , J of trapped modes 
where J is the greatest integer satisfying 

IT 

(40) (2J + 1) 6 < 2 . 

These higher discrete modes are algebraically complicated but retain longshore 

periodic behavior. Their crest amplitude is maximum at the shoreline but 

oscillates in the offshore direction with j nodes while decaying exponentially 

(see figure lib). To complete the set of normal modes, Hanson”' - '’ showed that 

there exists a continuous spectrum of waves with wave-number magnitude k 
2 

satisfying w > gk. These modes are even more complicated to describe 
algebraically but resemble simple deep-water wave trains far offshore. An 
example of one of these modes with crests parallel to the shoreline is shown 
in figure lib. We emphasize here that the continuous spectrum modes are not 
trapped like the discre'te modes. 
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Oceanographic interest in edge waves appears to have originated in 

Isaacs Zt aZ. 16 suggestion that these waves might be responsible for "surf beat" 

phenomena. Since then, edge waves have been indicted as a potential mechanism 

for numerous features of coastal dynamics such as beach cusps (Guza and Inman) 17 , 

crescentic bars (Bowen and Inman) 1 ^, and rip currents (Bowen and Inman) 1 ^. The 

lack of direct observational evidence of edge wave modes on beaches (i.e., wave 

crests pointing offshore!) naturally raises questions as to how and if these modes 

can be excited. The most important mechanism for extensive generation opportunities 

20 

of edge waves was provided by Guza and Davis who demonstrated theoretically 
that Stokes’ modes could be excited by simple wave trains normally incident onto 
a perfectly reflecting beach from deep water. Basically, the nonlinear coupling 
between the incident/refleeted wave of frequency W and background edge wave 
"noise" of subharmonic frequency u> = W/2 leads to a resonant interaction and 
growth of the edge wave noise. Eventually other nonlinear processes develop 
to limit growth; however, at steady state the edge wave amplitude exceeds that 

of the incident/reflected wave. The entire evolution of the edge wave has been 

. 21 
investigated theoretically by Guza and Bowen (hereafter referred to as G-B), 

22 23 

Minzoni and Whitman (hereafter referred to as M-W), and Rockliff 

Herein we present quantitative experimental measurements on edge wave 
excitation by normally incident wave trains. Some of the predictions of G-B and 
M-W are tested and necessary modifications for imperfectly reflecting beaches are 
presented. An outline of the presentation is as follows. In Sect. 2 we review 
the theoretical results with emphasis on the prediction which can be explicitly 
tested by the experiments. In Sect. 3 a brief description of the experimental 








facilities and data analysis techniques are discussed. A comparison of measured 
and predicted data appears in Sect. 4 followed by a summary of the major results in 
Sect. 5. More details of each aspect of this note may be found in Lin and Hammack.^ 
2. — Review of the Theories 

The first description of both edge wave excitation and subsequent evolution 

21 

was developed by Guza and Bowen using the shallow-water equations; hence, their 

results are limited to small beach angles (6 << 1). With clever but intuitive 

reasoning, G-B isolate several processes, analyze each separately, and combine 

23 

linearly to yield a complete evolution model. (More recently, Rockliff has 

reproduced some of the G-B results relying more formally on the governing 

22 

shallow-water equations.) Minzoni and Whitham use the full equations and provide 

the most formal (mathematically) description of the edge wave evolution process; 

their results remain valid for arbitrary beach slopes B £ ir/2. The greater 

formality of the M-W formulation permits more justifiable ordering of different 

processes. Hence, the outline below of edge wave evolution is based primarily 

on their formalism. Following the qualitative description, a quantitative listing 

of relevant parameters predicted by both M-W and G-B will be presented. 

The classical water-wave problem for the inviscid, irrotational, barotropic 

motions of an incompressible ocean in the wedge-shaped region of figure 11a may 

be formulated in terms of a complex velocity potential Initially, we have a 

linear wave train normally incident onto the beach from deep water; the incident 

wave amplitude far offshore is a^. We also take the incident wave to be perfectly 

reflected from the beach. (The reflection coefficient is R = a /a. =1 where a 

r l r 

is reflected wave amplitude.) Hence, the forcing for the onshore edge wave noise 
is a standing wave mode which is a member of the continuous spectrum of normal modes 









2 


discussed in Sect. 1. To simplify (somewhat) the standing wave description, 

it is further assumed that the beach slopes are a member of the denumerably 

infinite set 8 = ir/2N with N = 1,2, ... . The complex potential for the 

standing wave may be written as $ ' exp(iwt), and we note that ($ | ~ e << 1. 

In addition to the primary standing wave, we assume three-dimensional perturbations 

by Stokes edge wave with complex potential <j> e * exp(iwt) where | <J> | << e << 1 

initially. On a very short time scale (stage 1), nonlinear interactions between 

the primary wave and the edge wave noise are negligible. However, on a longer 

time scale (stage 2), quadratic interactions arise as a consequence of the nonlinear 

boundary condition at the free surface and become significant. Assuming the edge 

wave noise that satisfies the resonance condition w = W/2 to dominate, quadratic 

interactions <J> g <J> e * (where * denotes complex conjugate) give rise to terms of 

the form exp(iWt/2) and contribute to the growth of the subharmonic edge wave 

noise. This interaction between the linear standing wave and linear edge wave 

noise produces exponential growth of the form a ~ exp(yt) where a^ is the (real) 

edge wave amplitude at any crest location (i.e., a = a(x,y,t), and we take 

a(x,0,t) = A(x,t) as the shoreline amplitude) and y is the initial growth rate. 

The edge wave noise grows until a later time scale (stage 3) is reached where 

its finite amplitude leads to the development of processes which limit further 

growth. First, quadratic self-interactions of the edge wave with the form 

<|> 4> exp(iWt) contribute to the offshore standing wave; this is termed 

6 6 

radiation by G-B and corresponds to the fact that nonlinear (even second order) 
standing edge waves leak energy at frequency 2u> to deep water. Since 14> s | ~ c 
for the disturbed standing wave in stage 3, radiation feedback suggests that 
a steady state may be reached when |<J> e | ' e* 5 , i.e., the edge wave becomes larger 
locally than the forcing wave. Second, cubic self-interaction terms of the form 









t> <j> d> * exp(iWt/2) contribute further to edge wave nonlinearity by modifying 
e e e 

its dispersion relation. Hence, the system is retuned and the resonant frequency 
is shifted. If the forcing was perfectly resonant initially, it will now be off 
resonance and an effective reduction of the on-resonance growth rate y will occur. 
(Of course, if the initial forcing is slightly off-resonance, a similar modification 
in y is required.) When combined, all of these processes lead to an evolution 
for (inviscid) edge waves of the form 

dA 2 2 ** 3 

(41) — = (y - k ) A - y A 3 

where y is the initial on-resonance growth rate, k is a measure of the reduction 
in forcing efficiency due to nonlinear retuning, and y is a feedback coefficient 
due to both radiation and retuning. (Note that in (41) we tacitly assume that 
the edge wave phase remains constant through stage 3. Theoretically, this 
assumption is invalid; however, the experimental data support its applicability.) 

At the risk of misquoting G-B and M-W and possibly introducing numerical 
errors, we now list quantitative expressions distilled from their studies for 
the parameters appearing in (41). The initial (on-resonance) growth rate predicted 
for the edge waves is 

0.0424 a. W 3 

(42a) G-B; y = --—' 

g/T 8 

0.0426 a. W 3 

(42b) M-W y = -1-' 

g/S~ tan B 

where the M-W result is accurate to within 3% for B ir/4. These predictions 
are in agreement in the shallow-water limit (6 -*• 0) . Modification of the 
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initial growth rate by nonlinear retuning in steady state may be conveniently 

written to the same order of small parameter £ = ot k (where a is the 

0 0 0 

steady-state run-up amplitude of the edge wave) as 

.... „ ( 2 ) ( 0 ) 

(43) k = w - w 


with <i/ 2 ^ and representing the nonlinear and linear natural frequencies 

( 2 ) 

of the edge wave, respectively. The shallow-water limits for w given 
explicitly by G-B and inferred from the presentation of M-W are 


(44a) 

G-B: 

(2) 

0) 

(0) 

= 0) 

[i + 

0.055 e 2 ] 

0 

(44b) 

M-W: 

(2) 

U) 

(0) 

= . 0) 

[i + 

0.012 £ 2 ]; 
0 


= (g k B)* 5 in both cases. The different results of (44a) and (44b) do 
not appear reconciliable. Predictions for the feedback coefficient y in the 
shallow-water limit also differ at the same order of approximation according to 

W 5 

(45a) G-B: y = 0.001795 -- ■ 

g e 

(45b) N-W: u = 0.002323 

g s 


Finally, we note that the steady-state amplitude of the edge wave at the 
shoreline is found from (41) to be 

(46a) A 

o 


(Y 2 - 


2,4 

K ) 


L 


J 
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3. — Experimental Facilities and Procedures 

A set of five experiments are conducted to study the excitation of 
standing edge waves on a beach by waves normally incident from offshore. 
Experiments are performed in a laboratory basin 1.83 m wide, 4.0 m long 
and 30 cm deep as sketched in figure 12. The tank is equipped with a 
mechanical wave generator and a smooth beach whose slope can be varied. An 
actual beach width of b = 1.60 m was used in the experiments by incorporating 
internal sidewalls over the beach section as indicated in figure 12. This 
enables the offshore tank section of uniform depth to be surrounded by 
energy-absorbing material to minimize unwanted reflections and simulate 
open-ocean conditions. Extensive precautions were taken to insure that the 
shoreline and wavemaker were parallel and to prevent beach "pumping" during 
actual experiments. 

An offshore water depth of h = 25.4 cm was used in all five experiments 
while the beach slope was varied from 6 = 15° (N = 6) to 6 = 22.5° (N = 4) ? 
one experiment, run 4, was conducted at 6 = 20° which does not correspond to 
a slope of w/2N for any integer N. It is immediately apparent that the 
offshore uniform depth region does not conform to the (unrealistic) mathematical 
model where the water depth increases linearly offshor . v :ver, the offshore 
distance over the sloping beach section always exceeded three e-folding distances 
(> 95% decay) for the edge wave modes excited (assuming Stokes' modes). Based 
on Minzoni, 25 it is not expected that the edge wave dispersion relation is 
affected by the uniform depth offshore in any of the experiments. 

In all experiments wave amplitudes are measured at three locations; one 
gage onshore over the sloping beach and two gages offshore in the uniform 
depth region. The onshore gage is always positioned at a longshore location 
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corresponding to an edge wave antinode. The reality of edge waves very 
near the shoreline requires that they deviate from the profile predicted 
by (38); hence, the onshore gage is located approximately one e-folding 
distance away from the shoreline and always outside of the surf zone. 

The procedures adopted for each of the five experimental runs are 
as follows. First, the frequency w of the desired edge wave is calculated 
based on the beach width b = 1.6 m; the offshore wavemaker is then adjusted 
to a frequency W = 2w. The incident and reflected wave amplitudes (a^ and a^) 
of the partial standing wave systems are then determined from simultaneous 
measurements at the offshore gages in the absence of edge wave excitation. Edge 
wave excitation is suppressed by inserting a thin plate perpendicular to the 
shoreline which penetrates the surf zone; the plate introduces boundary 
conditions which destroy the resonance condition necessary for rapid edge wave 
growth. The thin plate is then removed permitting the evolution of the edge 
wave. Once the edge wave attains a steady state the wavemaker is stopped and 
the (unforced) decay of the edge wave over the sloping beach is measured. Hence, 
each experiment produces a time series of the water surface elevation onshore 
and offshore containing stage 1, stage 2, stage 3, and additional stage 4 where 
edge wave forcing is terminated and viscous damping forces dominate. 

4. — Comparison of Experiment and Theory 

Table 1 summarizes the measured data for the partial standing waves 
generated by the wavemaker/beach system in the absence of edge waves. The 
incident and reflected wave amplitudes are presented along with the respective 
reflection coefficients which range from 3 = 0.134 to R = 0.500. (Recall that 







theoretical results assume perfect reflection with R = 1.) It should be 

noted that wave breaking near the shoreline was observed in each experiment 

and was especially strong in run 1. 

A typical time series (run 2) taken onshore and offshore during edge 

wave excitation, evolution, and eventual damping is shown in figure 13. The 

onshore record clearly shows the effect of edge wave growth as it alternates 

between construction and destructive interference with successive crests of 

the incident/reflected wave system. (This behavior is a direct consequence of 

subharmonic excitation with W = 2a>.) A periodogram computed using the Fast 

Fourier Transform of the onshore wave record for run 2 is shown in figure 14. 

The first peak in the spectrum corresponds to the excited edge wave mode while 

the second peak corresponds to the standing wave generated from offshore. 

The equally spaced peaks at higher frequencies represent the superharmonics of 

both the edge wave and standing wave. The centering of wave content about 

well-separated and narrow bands with identifiable dominant frequencies (w, 2w, 3 oj, 

such as that exhibited in figure 14 permits the use of complex demodulation 

techniques to view the real time evolution of both the amplitude and phase of 

26 

each dominant wave component (see Bloomfield ). Using the periodogram, an 

initial guess, say w , for the dominant frequency in a band of interest is 

determined. The time series is then multiplied by exp(-iw t) to shift the 

a 

desired wave content to zero frequency. This demodulated signal is then 
low-pass filtered in time to yield the instantaneous phase and amplitude of the 
signal component with frequency cu . If the estimated frequency w is in error, 

g g 

a linear change in instantaneous phase with time will be observed; the slope of 
the linear change represents the error in In this manner the dominant 






frequency of each narrow band in the periodogram may be determined very 

accurately as well as its instantaneous amplitude and phase. (More details 

on the application of complex demodulation technique including a discussion 

24 

of the low-pass filter properties may be found in Lin and Hammack. ) The 

edge wave frequencies in each experiment were determined using this technique 

and are summarized in table 2. Figures 15 and 16 show the instantaneous 

amplitude and phase, respectively, of the edge wave harmonic at the onshore 

gage in run 2. Both parameters oscillate rather wildly during the initial 

time t < 25 secs; this behavior is a characteristic result when the signal-to-noise 

ratio is small. (Note that during this period this instantaneous amplitude 

of the edge wave is less than 0.07 cmi) The edge wave amplitude in figure 15 

then begins to grown in an exponential manner while the edge wave phase in 

figure 16 becomes constant. As time continues to increase, edge wave growth 

slows and a steady state is achieved for t > 70 secs; the gage-site amplitude 

of the edge wave at steady-state is a Q = 0.754 cm. The measured growth rates 

Y and steady-state amplitudes a are shown in table 2. The inferred shoreline 
in o 

amplitude A q based on the measured a^ at gage-site is easily calculated; the 
ratio of A^/a^ ranges from 3-6 in all experiments in agreement with the 
theoretical prediction that steady-state amplitudes of the edge wave should 
exceed those of the offshore standing wave. It is important to note that 
measured growth rates are significantly influenced by viscous and turbulent 
damping forces which must be clarified before a legitimate comparison with 
inviscid and irrotational theoretical models. The viscous damping rate is 
easily measured in the laboratory model simply by turning off the wavemaker 
and monitoring the damping of the edge wave harmonic with time. Typical 
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results are shown in figure 17 which corresponds to run 2. As expected, 

exponential decay occurs when the edge wave becomes sufficiently small (and 

linear theory is applicable). The measured decay rates 6 in the exponential 

stage for each experiment are presented in table 2. As noted by Guza and 
27 

Bowen, the damping of edge waves in the presence of breaking incident waves 

is likely to be dominated by turbulent exchange mechanisms rather than (laminar) 

viscosity. Viscous effects are enhanced by laboratory model scales and should 

play an important role; however, based on the results to follow, turbulent 

exchange mechanisms also appear to be important at laboratory scales. Ignoring 

the effects on laboratory damping of edge waves due to incident wave breaking, 

a measured "inviscid" estimate of the initial growth rate may be calculated 

by y ' = y + 6. Results for y ' are tabulated in table 2 along with the 

m ip m 

theoretical predictions for the growth rate according to 42a. It is evident 

that large discrepancies still exist; however, it should also be remembered that 

the theoretical predictions assume a perfectly reflecting beach. To examine the 

potential effect of the imperfect reflection in the experiments we have plotted 

the ratio of y ' /Y versus the reflection coefficient R in figure 18. Note that 

excellent correlation exists (data lie on 45° line) suggesting that the growth 

rates should be calculated using the reflected wave amplitude a^_ instead of the 

incident wave amplitude a^, i.e., we should replace y by yR- Theoretical estimates 

of yR are listed in table 2. It is important to emphasize that this behavior 

of the measured data clearly demonstrates that it is the reflected wave component 

of the offshore standing wave which drives the edge waves. It further supports 

17 

the observation (Guza and Inman ) that edge waves do not occur on dissipative 
beaches which produce little reflection. 
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In order to compute the theoretical amplitude of the edge waves in 
steady state, the results above suggest that we should modify the theoretical 
growth rate in (46) to the form 

2 ,*r 

R(y - o 

V 

The theoretical amplitude at a specific gage site then becomes 

(47) a = A cos k x exp(-k y cos B); 

oo m m 

results of computations for the predictions of both G-B and M-W are shown 
in table 2. Except for run 4, theoretical results exceed measured data by 
factors of 1.18 to 1.91. This behavior pattern is expected since the turbulent 
damping of the edge waves in the presence of breaking incident waves has been 
neglected. Also consistent with this hypothesis is the fact that agreement 
between experiment and prediction is best in run 5 where the least wave breaking 
was observed. Run 4 is anomolous as the measured amplitude exceeds theoretical 
prediction. (Interestingly, absolute agreement between measured and predicted 
data is best for run 4!) There is also the following evidence that run 4 is 
unusual re. ive to the other experiments. Note from table 1 that the offshore 
itiniina wave changes little between runs 3 and 4 while the steady-state amplitude 
• • wave in table 2 doubles! In other words, a major difference in 

■ ■ -:r- even thouah the forcing remains approximately constant. 


(46b) A = 

o 


: 


ween runs 3 and 4 is the slightly greater 
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beach slope (2°) in run 4 which no longer corresponds to one of the set 
it/ 2N required by the theoretical models. Based on this limited datum, it 
does appear that the predictions for steady-state response on beach with 
slopes satisfying 6 = t/2N cannot be extrapolated to intermediate slopes. 
(However, extrapolation for the initial growth rates does appear permissible 
with the modifications due to partial reflection described earlier.) 

5. — Conclusions 

Based on the experiments described herein, the following major conslusions 
may be stated regarding edge wave excitation on beac.nes by normally incident 
waves from offshore. 

i) The theoretical growth rates calculated for perfectly reflecting beaches 
may be generalized to beaches with imperfect reflection simply by 
reducing the growth rates in direct proportion to the reflection 
coefficient of the beach. This result suggests that the reflected wave 
is indeed the driving force for edge wave response and absolutely 
necessary in order to excite these modes. 

ii) Growth rates modified for partial reflection effects may be extrapolated 
to beaches with slopes not satisfying S = tt/2N. 

iii) No change in the edge wave phase occurs once the exponential growth 
stage is encountered. Hence, the evolution equation for the real 
amplitude of the edge wave is applicable. 

iv) There is limited evidence that the steady-state amplitudes of edge 
waves on beaches with 8 7* ff/2N are significantly larger than those 
for nearby beach slopes with 8 = t/2N. 


i 





v) Damping effects on edge waves are significantly influenced by turbulent 
exchange mechanisms resulting from breaking incident waves, even on the 
laboratory scale. 
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FIGURE CAPTIONS 


Figure 1. Definition sketch of the fluid domain. 

Figure 2. Schematic drawing of the wave generator. 

2 

Figure 3. Surface wave evolution: h = 5 cm, b = 61 cm, V = 30.5 cm , 

N = 3. ' ■■ , measured profiles; # , soliton profiles 

computed using (7) . 

2 

Figure 4. Surface wave evolution: h = 5 cm, b = 30.5 cm, V = 30.5 cm , 

N = 3. «— , measured profiles; # , soliton profiles 

computed using (7). 

2 

Figure 5. Surface wave evolution: h = 5 cm, b = 61 cm, V = 30.5 cm , N = 0. 

2 

Figure 6. Surface wave evolution: h = 5 cm, b = 30.5 cm, V = -30.5 cm , N = 1. 

2 

Figure 7. Surface wave evolution: h * 5 cm, b = 30.5 cm , V = 0, N = 2. 

— h i ■ ■ , measured profiles; • , soliton profiles computed using (7). 

Figure 8. Internal wave evolution: h^ = 45 cm, h 2 =5 cm, 
b = 61 cm, A = 0.05 

Figure 9. Comparison of leading wave profiles with theoretical 
(internal) KdV soliton. 

Figure 10. Evolution of a narrow-banded wave packet. - » , measured 

profiles; # , soliton profile computed using (12). 

Figure 11a) Definition sketch of fluid domain and Stoke's edge-wave mode. 

b) Sample offshore profiles of higher edge-wave modes: —— , 

discrete spectrum (trapped) mode; ——, continuous spectrum 
(untrapped) mode. 







I 


- 40 - 


Figure 12. 

Figure 13. 
Figure 14. 
Figure 15. 

Figure 16. 

Figure 17. 

Figure 18. 


Schematic drawing of wave basin — plan view. 

Time-series of water surface elevations onshore and offshore. 

Periodogram of onshore time series. 

Amplitude evolution of subharmonic edge wave for run 2. 

Phase evolution of subharmonic edge wave in run 2. 

Viscous damping of subharroonic edge wave in run 2. 

Correlation of measured and theoretical "inviscid" growth 
rates with beach reflection coefficient. 
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